Beyond the Quantum Adiabatic Approximation: Adiabatic Perturbation Theory 
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We introduce a perturbative approacli to solving the time dependent Schrodinger equation, named 
adiabatic perturbation theory (APT), whose zeroth order term is the quantum adiabatic approxi- 
mation. The small parameter in the power series expansion of the time- dependent wave function is 
the inverse of the time it takes to drive the system's Hamiltonian from the initial to its final form. 
We review other standard perturbative and non-perturbative ways of going beyond the adiabatic 
approximation, extending and finding exact relations among them, and also compare the efficiency 
of those methods against the APT. Most importantly, we determine APT corrections to the Berry 
phase by use of the Aharonov-Anandan geometric phase. We then solve several time dependent 
problems allowing us to illustrate that the APT is the only perturbative method that gives the right 
corrections to the adiabatic approximation. Finally, we propose an experiment to measure the APT 
corrections to the Berry phase and show, for a particular spin-1/2 problem, that to first order in 
APT the geometric phase should be two and a half times the (adiabatic) Berry phase. 

PACS numbers: 31.15.xp, 03.65.Vf 

Keywords: Perturbation theory; Phases: geometric; dynamic or topological 



I. INTRODUCTION 



Aside from interpretation. Quantum Mechanics (QM) 
is undoubtedly one of the most successful and use- 
ful theories of modern Physics. Its practical impor- 
tance is evidenced at microscopic and nano scales where 
Schrodinger's Equation (SE) dictates the evolution of the 
system's state, i.e., its wave function, from which all the 
properties of the system can be calculated and confronted 
against experimental data. However, SE can only be ex- 
actly solved for a few problems. Indeed, there are many 
reasons that make the solution of such a differential equa- 
tion a difficult task, such as the large number of degrees 
of freedom associated with the system one wants to study. 
Another reason, the one we want to address in this pa- 
per, is related to an important property of the system's 
Hamiltonian: its time dependence. 

For time independent Hamiltonians the solution to SE 
can be cast as an eigenvalue/eigenvector problem. This 
allows us to solve SE in many cases exactly, in particular 
when we deal with systems described by finite dimen- 
sional Hilbert spaces. For time dependent Hamiltonians, 
on the other hand, things arc more mathematically in- 
volved. Even for a two-level system (a qubit) we do not, 
in general, obtain a closed-form solution given an arbi- 
trary time dependent Hamiltonian, although a general 
statement can be made for slowly varying Hamiltonians. 
If a system's Hamiltonian H changes slowly during the 
course of time, say from t = to t = T, and the system 
is prepared in an cigcnstate of H at t = 0, it will rc- 
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main in the instantaneous (snapshot) cigcnstate of H(t) 
during the interval t G [0,r]. This is the content of the 
well-known adiabatic theorem [l|. 

But what happens if H(i) is not slowly enough varied? 
For how long can we still consider the system to be in a 
snapshot cigcnstate of H(i), i.e., for how long the adia- 
batic approximation is reliable? What are the corrections 
to the adiabatic approximation? One of our goals in this 
manuscript is to provide practical and useful answers to 
these questions. We introduce a perturbative expansion 
about the adiabatic approximation, named adiabatic per- 
turbation theory (APT), using the quantity v = \/T as 
our small parameter. This power series expansion in v 
is subsequently used to calculate corrections to the adia- 
batic approximation for several time dependent two-level 
systems. It is worth noting that answers to previous ques- 
tions can also be seen, under certain provisos, as a way of 
solving perturbatively any time dependent problem. We 
should stress that the APT is not related to the time- 
ordered Dyson series method since the latter is not a 
perturbative expansion about the adiabatic approxima- 
tion, in terms of the small parameter v. Rather, it is an 
iterative way of getting the unitary operator governing 
the evolution of a system, in terms of a small perturba- 
tive potential in the Hamiltonian. 

Another goal is to present an exhaustive comparison of 
all the approximation methods developed so far to solving 
SE. In particular, we show the exact equivalence between 
Garrison's multi-variable expansion method [4] (which 
solves an extended set of partial differential equations) 
and APT. However, it is important to stress that the 
APT, being an algebraic method, is straightforward to 
use while Garrison's approach is very hard to extend be- 
yond first order. We also provide an extension to Berry's 
iterative method [3| where, contrary to the original ap- 



proach, we keep all terms of the new Hamiltonian ob- 
tained after each iteration. We then discuss the possibil- 
ity to choose other types of iteration (unitary transforma- 
tions) to potentially do better than Berry's prescription. 

Furthermore, it is known that if the conditions of the 
adiabatic theorem are satisfied and H(r) ~ H(0), it fol- 
lows that the state |^(T)) describing the system at i = T 
is given by |^(T)) = e''^(^)|1'(0)), where |^(0)) is the 
initial state and (j){T) is a phase that can be split into 
dynamical and geometrical parts [4| • This raises another 
question we address here and which is not independent 
from the ones above: what are the corrections to the 
Berry phase [j| as the system deviates from the adia- 
batic approximation? To provide an answer we make 
use of the Aharonov-Anandan ( AA) geometric phase Q , 
which is a natural extension of the Berry phase having a 
geometric meaning whenever the initial state returns to 
itself, even for a non-adiabatic evolution. We thus com- 
pute the AA phase for the corrections to the adiabatic 
approximation which, therefore, possess the geometrical 
and gauge invariancc properties of any A A phase. We 
then show, for a particular spin-1/2 example, that when- 
ever H(T) = H(0) and the evolving state corrected up to 
first order returns to itself (up to a phase) at i = T, we 
obtain a geometric phase that is two and a half Berry's 
phase value. 

In order to provide a clear and complete analysis of the 
questions raised above we structure our paper as follows. 
(See Fig. [T]for a structural flowchart of the paper.) In 
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FIG. 1: Different approximation methods to solving the time- 
dependent Schodinger equation. APT: Adiabatic pertur- 
bation theory (Garrison, Ponce, this paper); IRBM: Itera- 
tive rotating-basis method (Kato, Garrido, Nenciu, Berry); 
TDPT: Time-dependent perturbation theory (Dirac); SA: 
Sudden approximation (Messiah); AA: Adiabatic approxima- 
tion (Born and Fock). 

Sec. lUwe review the adiabatic approximation, highlight- 
ing the conditions that the snapshot eigenvectors and 



eigenvalues of H(i) must satisfy for this approximation 
to be valid. In Sec. IIIII we review many strategies that 
may be employed to find corrections to the adiabatic ap- 
proximation as well as to the Berry phase. As shown 
later, those methods are unsatisfactory since either they 
do not furnish all the terms that correct the geometrical 
phase and the adiabatic approximation or they cannot 
be seen as a perturbation in terms of the small param- 
eter V ~ IjT. In Sec. IIVI we present our perturbation 
method, i.e. APT, in its full generality and provide ex- 
plicit corrections to the adiabatic approximation up to 
second order. In Sec. |V]we deal with corrections to the 
geometric phase using the previous method, presenting 
its first order correction. In Sec. IVII we compare all 
other methods with the APT, emphasizing the main dif- 
ferences among them. In Sec. IVIII we review the exact 
and analytical solution of a time dependent problem and 
expand it in terms of the small parameter v. Then we 
show that our perturbative method is the only one that 
gives all the terms obtained from the expansion of the ex- 
act solution. We also propose an experiment where APT 
corrections to the Berry can be measured. In Sec. IVIIII 
we solve numerically three other time dependent prob- 
lems and compare them with our perturbative method. 
Finally, in Sec. IIXI we provide our concluding remarks. 



II. THE ADIABATIC APPROXIMATION 

Let us start rewriting the time dependent SE in terms 
of the rescaled time s — vt, where T = 1/v is the relevant 
time scale of our Hamiltonian H(i). We then formally 
solve the SE, emphasizing the assumptions imposed on 
the spectrum of H(i), and show the conditions the in- 
stantaneous (snapshot) eigenvectors of H(i) must satisfy 
for the adiabatic approximation to be valid. 

The time dependent SE is written as 



in-\^{t))^ii{t)\^{t)), 



(1) 



where \'^{t)) is the state describing our system at time 
t. Since we want to work with the rescaled time s and 



^ = w ^ it results 



i;i«-|*(s))=H(s)|vI/(s)). 



(2) 



Building on the knowledge that the adiabatic phase can 
be split into a geometrical (7) and a dynamical (uj) part 
[j| we may write down the solution |^(s)) as 

n=0 

in which 6„(s) are time dependent coefficients to be de- 
termined later on. The sum over n includes all snapshot 
eigenvectors of H(s), 



Uis)\nis))=E^{s)\n{s)), 



(4) 



with eigenvalue -E'„(s) (n = represents its ground state 
(GS)). The Berry phase associated to the eigenvector 
|n(s)) is 



7„(s)=i/ {n{s')\—n{s'))ds' = i I Af„„(s')ds', (5) 



while 



w„(.s) 



En{s')ds' = VUJn{t) 



(6) 



defines its dynamical phase. Let us start assuming that 
H(s) has a non-degenerate spectrum during the whole 
evolution. Note that the initial (s = 0) conditions on 
|^(s)) are encoded in 6„(0). Therefore, if the initial 
state is |0(0)) we will have 6„(0) = (5„0j where Stj is the 
Kronecker delta. In this case, as we will see below, the 
spectrum needs to satisfy the less restrictive condition 
i?o(s) 7^ En{s), Vs € [0,T], n 7^ 0, for our perturbation 
method to work. In other words, our method will work 
whenever one starts the evolution at the GS and there 
is no level crossing between Eo{s) and any other En{s) 
(even though the excited state part of the spectrum may 
display level crossings). Similar type of conditions can 
be shown to apply to states living in subspaces spectrally 
separated from the rest. 

Replacing Eq. ^ into (O) using Eq. ^ and left mul- 
tiplying it by {m{s)\ leads to 

b„{s) + ^ e->'""(^)e"''""(^)M„,„(s)6,„(s) = 0, (7) 

rn—O 

where the dot means 4- and the indices m ^-^ n were 

as 



exchanged. Here u;„ 

7™(s) ^7n{s)., and 



i(s) = U;m{s) - W„(s), 7m„(s) 



Mnrnis) = (n{s)\m{s)) . 



(8) 



So far no approximation was invoked and in principle 
the time dependence can be found by solving the system 
of coupled differential equations given in ([7]). General 
numerical methods to solve such equations will face the 
computational difficulty of integrating highly oscillatory 
terms such as e~"'^""'(*)e'''''""^*\ making the approach 
numerically unstable. Later on we show that our pertur- 
bative method gets rid of this problem. 

The adiabatic approximation consists in neglecting the 
coupling terms (O, i.e., setting A/„,„(s) = 0, 

bn{s) = bn{0) — > adiabatic approximation. (9) 

Replacing Eq. ([9]) into ([3]) we obtain, 

l*'°'(s)> = E e'^"(^)e->"(-^)6„(0)|n(s)), (10) 

where we used [^^"^(s)) instead of |^(s)) since the adi- 
abatic approximation will be the zeroth order term in 



the pcrturbative method developed later. In the case the 
system starts at the GS, 



*(°)(s)) 



ei7o(s)g- 



^"«('*)|0(s)). 



(11) 



For the sake of completeness, let us analyze some gen- 
eral properties of Mnm{s)- Since the eigenvectors of H(s) 
are orthonormal we have (n(s)|m(s)) = Snm- Taking the 
derivative with respect to s we get Af„,„(s)-|-Mj^„(s) = 0, 
which implies that A'/„„(s) is a purely imaginary number, 
as it should be since 7n(s) is real. When n ^ m, by tak- 
ing the derivative of Eq. Q with respect to s and left 
multiplying by {m{s)\ one gets 



M„,„(s) = (7i(s)|H(s)|m(s))/A„„(s), 



(12) 



where Am„(s) = Em{s) — £'„(s). This last expression 
indicates that the adiabaticity condition is related to the 
existence of a gap. A spectrum of discussions on the 
validity of the adiabatic approximation can be found in 
Refs. 
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III. CORRECTIONS TO THE ADIABATIC 
APPROXIMATION 

We can classify all the strategies to find corrections to 
the adiabatic approximation into two groups. The first 
one includes those methods that perform a series expan- 
sion of the wave function in terms of the small parame- 
ter V = 1/T <^ 1, with T representing the time scale for 
adiabaticity. In this group we include the pioneering ap- 
proach of Garrison ^ and the seminal work of Ponce et 
al. [ll|. The second group includes those methods that 



intend to approximate the solution to the time depen- 
dent SE without relying on a formal series expansion of 
the wave function 0, [l3, Ej 13 ^^^ using the adiabiatic 
approximation as their zeroth-order step. In this section 
we review two methods belonging to the first group and 
one to the second, called adiabatic iteration by Berry [3|. 
We then comment on a possible extension of the latter. 



A. Examples of the first group 

We first show how to manipulate Eq. ([7]) in order to 
get a series expansion in terms of the small parameter d, 
which we call the standard (textbook) approach. We then 
discuss the multi-variable expansion method of Garrison 
0, who also dubbed it APT. 



1. The standard approach. 
One can formally integrate Eq. ^ to obtain 
bn{s) = bn{0) -Y. f ds'e->™"(^')B,„„(s'), (13) 

m=0 -^0 



where 



5„m(s) =e'^'""(^)M„,„(s)fe,„(s). (14) 

The integral inside the sum in Eq. (jl3|) can be written as 

'ds'i?„„(s')e^^°''i^"^'""(^"\ (15) 



in which Cmn{s) = —iAmn{s)/fi. Our goal here is to 
expand / in powers of v. This can be done by using the 
mathematical identity 



B™(.)e^/o d.'c„„(.') ^ d / ^™«(£)^i/=d.'c_(.')') _ d fB^\ ^i;„m.'c_(s'). (^g) 

as \ Cjnn{s) J as \Cmn\s) J 

Replacing Eq. pB]) into ([15]) we arrive at 

U™nW C™„(0)y io ds' lvC„4s')r ■ ^ ' 

One can apply the identity (|16p again to the integrand of the last term by substituting Bmn{s) for v-^ ( ^-j"'";^^ ] , 



\^mn(s) 



Smn(O) 
L-mnCJ) 



1 



t^mn(Sj ds V 0^)1(5 j 
I 



gmn(g) \ 1 f„-d.'C„„(.') 



C(t'^), (18) 



with the symbol 0{v^) standing for the term 



replacing hm{s) —> &m(0) at the right-hand side of ([T 
bnis) = 6„(0) - ihv 



■■2 ' ds'- 



ds' \Cmn{s') ds' \Cr„n{s') 



i/o ds"C,„„(s") 



E 



^--^^As)^h^,ds)Ml^b^^Q) 



^mn \S ) 



One can similarly continue the iteration to obtain higher ^ ' 

order terms but the first two are already enough for our Finally, substituting Eq. (gOl) into Q we get the (un- 

purposes. We should note that, strictly speaking, the normalized; normahzation introduces higher order cor- 

procedure just described is not a genuine power series rections in v) state that corrects the adiabatic approxi- 



mation up to first order via the standard approach, 

|*(s)) = |*W(s)) + vl^^^s)) + 0{v^), (21) 



expansion in terms of the small parameter v. This is be- 
cause to all orders we have a phase contribution {Cmn{s) 
is purely imaginary) of the form e~ Jo ^'^ Cmn{s )_ xhis 
term is related to the dynamical phase of our system and where |^(°)(s)) is given by Eq. pO|) and 
together with the Berry phase will play an important role 
in the APT developed in Sec. IVl 



Using Eq. ([T8|) in (fT3|) and keeping terms up to first or- 
der in V we obtain after substituting the values of Bmn{s) 
and Cjnnis) 



bn{s) = &„(0) - ihv 



E 

m— 
m^7i 



|*(i)(s)) = ihY^ e-^ 

n,?7i— 

-i^E 

n,m— 






As)^i 



6™(0)|n(s)) 



|M„„(0), 



g i^^i^)e''r^(^ yr"'):; buo)Hs)), 



A„™(0) 



.is)^n„Ms)Ml}I!^b^(^s) 



^mn\^ ) 



(22) 



(19) 



Note that we have to solve this equation itcratively keep- 
ing terms up to first order in v. This is equivalent to 



with A,„,i(s) = — A„,„(s). If the system is at the GS at 
s ~ Q, 6„(0) = 5nO; and Eq. (P^ reduces to 

|*(i)(s)) = iftye->°(^)e'^°(^)^^^|n(s)) 

- iftE^~"""^^^«"^"^^^X^I-(^))'(23) 



which displays no hncar in v correction to the |0(s)) com- 
ponent (the sum starts at n = 1). As shown in Sec. IIV[ 
there is a missing term correcting the coefficient multi- 
plying the GS that naturally appears in the APT. Also, 
|^(^)(0)) = 0, as we would expect since we must recover 
the initial state |*(°)(0)) at s = 0. 



To solve Eq. ((24|) we write the wave function as follows 

|vI/(s))=5^e->"c„Ks)Ks)), (25) 

n=0 

where uj represents all the variables w„ and 



2. Multi-variable expansion method 

To obtain a time dependent multi-variable SE we con- 
sider the quantities u>nis) as independent variables, i.e. 
^n{s) — *■ w„ 1^. They are called fast variables in con- 
trast to the rescaled time s, which is the slow variable. 
In this language the differential operator v^ is replaced 
by vds 



-Du,, where 



D„ 



E, 



Eir"^-"' 



n=0 



and the modified SE is written as, 

ihivds + D^.)\^{s)) ^il{s)\'^{s)). 



(24) 



c„(c.,s) = ^z;''c(f)(a;,s). 

p=0 



(26) 



Note that c„(tj, s) is written as a power series in v and our 
goal is to obtain Cn {w, s) to all orders. Using Eq. ((26|) 
we can rewrite (l25l) as 



l*W) =EE«'«'"""^»'^('^'^)I"(^))- (27) 



Substituting Eq. ^ in the modified SE (Eq. ([24])), car- 
rying out the derivatives, and taking the scalar product 
with {m{s)\ we get 



p=0 \ n=Q ) p=0 



Noting that the last term of the previous equality can be 
written as 

oo 

Y^ yPe-^^'- D^cI^\lu, s) = e-^'-^D^c^°\uj, s) 

p=0 

we can rewrite Eq. (P5]l in the following form 

oc 

p=0 

+ ^e-^--M„„(s)c(^)(^,s)J 

m=0 / 

+e-"'""D„ci°)('^,s) = 0, (29) 

where we have exchanged n <-!■ to. A sufficient condition 
for the validity of Eq. ^^ is obtained when we set 



D^c^J;'>{io,s) = 0, 



(30) 



and 



+ ^c->-M„™(s)c(Pnw,s)=0. (31) 



r 



Hence, we can calculate the coefficients c„ (w, s) by solv- 
ing the partial differential Eqs. ([50]) and ([5T|) . Note that 
to seek for the solution of order p we need to have the pre- 
vious, p—1, order solution. Furthermore, as we increase 
the order, the partial differential equations become more 
cumbersome constituting a practical limitation of this 
method. The APT developed in Sec. IIVI on the other 
hand, does not rely on any differential equations what- 
soever. All corrections to the adiabatic approximation 
of order p are obtained via algebraic recursive relations 
that involve coefficients of order p—1. This will allow 
us to derive in a relative straightforward manner explicit 
expressions up to second order in the small parameter v. 

In what follows we derive explicit expressions for 
c„ (w, s) and c„ (w, s). To zeroth-order Eq. ((30)) tells us 
that Cn {i^,s) does not depend on the variables uj, i.e., 
c„ ((jj,s) = c„'{s). Moreover, since at s = we have 
the initial condition |*(0)) = J2n=Q^n{0)\n{0)) then it 

immediately follows that c„ (0) = 6^(0) and 



„(p) 



(0,0) = 0, 



P^O. 



(32) 



m=0 



To have the adiabatic approximation as the zeroth or- 
der term in the power scries solution we must have (cf. 



Eq. §1 with dill)) 

cio)(s)=e"'"(^)6„(0), (33) 

which according to Eq. (pij) leads to 

D„c«(c^, s) + dj°\s) + M„„(s)ci")(s) 

+ ^ e->"-M„,„(s)c(^) (s) = 0. (34) 






But Eq. ([33]) together with jS]) imply that dsc'i?\s) 
M„n{s)c^n\s) = 0. Thus, Eq. dM]) becomes 



i^^cii'(w,s) + Y, e-^"-e'T'"(^)Af„„(s)&„(0) = 0, 



ni—0 



and wc now want to solve this equation. 
Following Garrison [2| we write 



(35) 



c(f)(a;,s) = c(f)(s)+d(f)(c.,s), (36) 

with the assumption that (average over uj) 

{d^P\u,s))^ ^ {D^d^P\uj,s))^ = 0. (37) 

In other words, we have separated out the u and s de- 
pendence of Cn into two contributions; the first depends 
only on s, and is called the average term; the second one 
depends on both w and s, but with the additional con- 
dition that its average over the fast variables w is zero. 

Thus, (ci^'(s))^ = c'^\s). Substituting Eq. ^ into 
(ESI) we get 



D^d\l\oj,s) + Y. e-""'""c'^-(^)M„„(s)6™(0) = 0, 

(38) 



m— 



and solving for dn we obtain 
d«(c.,.) = ift ^ c->'""c'-'"(^)^=i4^™(0)- (39) 

771 7^ n 

Note that dn {lo, s) +a{s), with a{s) independent of the 
variables w, is also a solution of Eq. ([35]). However, since 

we imposed that {dn [uJ, s))^ = 0, the only possible value 
for a{s) is zero. 

If the initial state is |0(0)) (6„(0) = (5„o) one gets 

d,W(a;, .) = ife-t-o„ei7o(.)^Zio(f) (^ _ ^^^^ (40) 

and since (d„ (tj, 5))^^ = and the only dependence on 
UJ in Eq. (|39|) is in c""'^"" we get 



We are now able to determine the average term Cn (s). 
Inserting Eq. ([55)) into ([3T|) we get for p = 1, 

Z?^di2)(w, s) + 934i)(u., s)+a,c,W(s) + M„„(s)4i'(w, s) 
+M„„(s)cii)(s) +Y. e~""'""M„™(s)d,W(c^, s) 

7717^71 

+ ^e-^"'""M„,„(s)cW(s) = 0, 

m=0 

where we have used that D^^Cn \s) ~ 0. Averaging over 
w, and noticing that {D^jdn (w, s))ui = (dL {i^, s))^^ = 0, 
(a,di^^(a;,s))^ = ds{di!\uj,s))^ = 0, and using Eq. ([SD 
we obtain 

dsci^Hs) + Mnn{s)c^n\s) 

+ Y Mn™(s)(c-""'""4'n^,s))„ = 0.(42) 

m=0 

We can recast the average (using Eq. p9|) ) as 

(e-t-.rfW(^,,))^ = i;,^^!!!^e^7.(.)(e-t-..)^5fc(0) 
;^QA„fe(s) 

= i?,^!!^e'^"(^)6„(0), (43) 

in which wc have used that i^mn + w/cm = '^fcn- Equa- 
tion ([32]) plus M„„(s) = -M;„„(s) imply that Eq 
can be written as 



ds 



+ p{s)cW(s)^qis), 



where 



(44) 



(45) 



p{s) = A/n,i(s), 

lis) ^ift^ l^-(;)l% i7.(-)5„(0), (46) 

and whose well known general solution is 
1 



cL'Hs) = 



^,{s')q{s')ds' + cl'\0)\, (47) 






It is interesting to note that the integrating factor ^{s) 
is related to the Berry phase 7n(s). Inserting Eqs. ([^5]) 
and ([46]) into ([48]) we get 



5W(.) = i;.e-(^) rd.'^^^^4^5„(0) 



m=0 



(e->-)^ = 5„ 



(41) 



,i7.(s)F(l) 



-ci!'m- 



(48) 



We can now write down the expression for c„ given d„ 
(Eq. §^) and cl^'> (Eq. (HH])), 



+iW^"(') 



^w„„„i7„(s) 



Mm«(s) 



^„(0) 



Jo n 






m^n 



gi7.(^)c(l)(0). 



(49) 



To determine c„ ^ (0) we use Eq. ([5^ , which guarantees 
that the adiabatic approximation is obtained as zeroth 
order, 



cW(0) = -ift^ 



m=0 



MnmjO) 
A„,„(0) 



bmiO). 



(50) 



Finally, expressing |^(s)) as given in Eq. ([2T|) and us- 
ing Eqs. ((271), ([49]), and ([SOl) we get for the first order 
correction to the adiabatic approximation, 

|vl/(i)(s)) = ihJ2 e->"(^)e'^"(^)J™„(s)6„(0)|n(s)) 



n,m— 
n,rri— 






.(s)pi 



&,n(0)|n(s)) 



- iS^ c-t-(^)e^-"(^)^J^fe,„(0)|n(.)), 



n,?n— 



in which 



'Jrtxn \ ^ } 



l\\2 



^ , ,|Mmn(s')| 

— A } — T 

A„i„(s ) 



(51) 



(52) 



Note that now we are writing again explicitly the depen- 
dence of Wn on time, i.e., uJn -^ LOn(s). For completeness, 
we write down the first order correction when we start at 
the GS (6„(0) = Sno) 

|^(i)(s)) = i;i^e-^"«(^)e'^«(^)j„o(s)|0(s)) 

n=l 

+ iftye-^"«(^)e'^«(^)^^^|n(s)) 
^^ A„o(s) 



n=l 



:M„o(0), 



ift^e-t-(^)e'^"(^):p^|n(.)). 



n=l 



A„o(0) 



(53) 



where we have replaced m —^ n in the first sum. 

Comparing Eqs. ^5^ and ^5^ with Eqs. (gj) and 
we immediately see that now we have a new extra term 
for the first order correction, the one proportional to 



Jmn{s). We would like to remark, though, that in Gar- 
rison's original work [2] he only obtained the first line in 
Eq. ([53|) . and thus our presentation constitutes an elab- 
oration on his general idea. Going beyond first order in 
V within Garrison's approach is an extraordinary tour de 
force. Fortunately, we will sec in Sec. IIVI that not only 
the extra term appears in our APT but, moreover, it is 
quite easy to obtain higher order corrections. Indeed, we 
will prove the mathematical equivalence between the two 
methods. 



B. Example of the second group 

The iterative method proposed by Berry [3[ consists 
of successive unitary operations that hopefully rotate the 
original basis or axes (the eigenvectors of the original 
Hamiltonian) closer and closer to the evolving state. In 
the most optimistic scenario a finite number of rotations 
would bring us to a moving frame in which the Hamil- 
tonian, as seen from this new frame, becomes time inde- 
pendent (this is the case in the simple single spin problem 
of Ref . [21| ) . Then we can solve the transformed Hamil- 
tonian using well developed time independent techniques 
and, by reversing the transformations, we would have the 
answer to the original problem. 

Berry [3| was only interested in corrections to the geo- 
metric phase that can be obtained by such a procedure. 
He showed that this strategy leads to successive correc- 
tions to the Berry phase although only in an asymptotic 
sense, i.e., after, let us say, the k-th rotation, the next 
following terms cannot improve the result achieved up to 
this iteration; rather, they spoil any possible useful cor- 
rection. In Ref. [3[ it was also shown, and we will review 
it here, that this iterative process is not an expansion in 
the small parameter v since every iteration contains v to 
infinite orders. We should also note that, as stated in 
Ref. lj|. Berry's iterative method is equivalent to the 
onesofRefs. [ll,[Il,[ll. 



In what follows we will extend Berry's approach to in- 
clude corrections to the wave functions. For the ease of 
notation, and since we will be dealing with successive 
iterations, we will denote the original Hamiltonian, its 
eigenvalues, and eigenvectors as H(°)(s), £"„ (s), and 
\n^^'{s)), respectively; after j iterations we will have 

H(-')(s), En (s), and \n^^\s)). Also, as in previous sec- 
tions, the initial state is written as [^'■"•'(O)). 

The main idea behind Berry's approach lies in the re- 
alization that the unitary operator Uo(s) (Uo(s)Uj(.s) = 
Uj(s)Uo(s) = 1) that gives the snapshot eigenvector of 
H(o)(s),i.c., 



|n(o)(s))=Uo(s)|n(o)(0)), 
can be used to construct the state 

|vI/W(,s))=Ut(.)|vI/(0)(5)), 



(54) 



(55) 



whose time evolution is determined to be 

i;iw|*(i)(s)> =H(i)(s)|*(i)(s)>, (56) 

with 

H(i)(s) = Ut(,s)HW(,s)Uo(s) - ihvljlis)iJois). (57) 

Repeating the previous argument with a new unitary op- 
erator Ui(s), which gives the snapshot eigenvectors of 

H(i)(s), 



|nW(s))=Ui(s)|n(")(0)), 



(58) 



ahows us to generate a new state |4'(^)(s)), and by iter- 
ating this procedure j times we obtain 

= U]_,{s)lj]_,{s) ■ . .Ul(.)Ut(.)|vI/(")(,s)), 
that satisfies the SE 

ihv\¥^^s)) = H(^')(s)|*(^')(s)), (59) 

with \n^3\s)) =Uj(s)|nW(0)) and 

HW(s) = Ut_,(s)H(^--i)(s)U,-i(.) 

-invU]_,(s)U,_i(s). (60) 

Using that M^n^\s) = {m^^-'^'>{s)\n^^~^\s)), the ma- 
trix elements of H^^^s) are 

(mW(0)|H(^)(.)|n(")(0)) = Ei=~'Hs)Sn,n~ihvM^;^'H-^)- 

(61) 
Loosely speaking, |\I'''')(s)) can be seen as the state ob- 
tained after cancelling or freezing (Uj_ j (s)) the time evo- 
lution of the snapshot eigenvectors of H^-'^^'(s), i.e., we 
are always trying to suppress the time dependence of the 
new Hamiltonian H'-''(s). Before we move on we should 
remark that Uj(s) is not the usual unitary operator Z-/j(.s) 
that evolves an arbitrary state l^'-'-'(O)) into the state 
|*(j)(s)), i.e., |*(^Hs)> =^j(s)l*'''HO)). 

Let us now explicitly show how to determine the state 
|\E'^-''(s)) [la]. For this purpose we write it as 

\^U)(s)) = 5]e"'^^"''(^)e->i^^"''(^)&«(s)|nW(0)), 

n=0 

(62) 
in which 7„ '(s) = i J^ Mnn (s')ds' is Berry's phase 
for the snapshot eigenvector \n^^~^^{s)), with dynami- 
cal phase uji^~^'{s) = ^ /p" E'n ~^ (s')ds'. Note that as 
opposed to Eq. ([3]), the eigenbasis used in (|S^ is not 
changing over time, i.e., instead of the snapshot eigen- 
vectors \n^^\s)) we now have |n(°^(0)). But as before, 

our goal is to find the equations satisfied by 5„ (s) which 
are obtained after inserting Eq. (|62p into ([59|) : 



m=0 



where UJmn (s) = USfi (s) - UJn (s) and Jmn (s) = 

7m"^^(s) - ji^~^\s). We see that Eq. §^ is formally 
identical to Eq. ([T]), which means that any technique de- 
veloped to solve (O can be employed to solve ([55]) : in 
particular the APT of Sec. IIVI Moreover, this formal 
similarity between these two equations evidences that 
Berry's iterative procedure is not a pcrturbative expan- 
sion about the small parameter v. Actually, as already 
anticipated, after each iteration we still have (in general) 
terms involving v to all orders. 

In closing, let us indicate a way to, in principle, ex- 
tend Berry's iterative approach. One can easily check 
that unitary iterations not constrained by the relations 
\n^^>[s)) = Uj(s)|n'°^(0)) lead to the same formal set of 
equations previously derived. Nonetheless, for a given 
number of iterations, the optimal choice of unitarics ap- 
proximating the real time evolution is a difficult problem 
related to the complexity of efficiently approximating an 
arbitrary unitary operator in a quantum circuit. 



IV. ADIABATIC PERTURBATION THEORY 

The reasons for introducing an APT are three fold. 
First, APT is a method that allows straightforward eval- 
uation of corrections to the geometrical phase (Berry 
phase). Such corrections are presented as a power se- 
ries in terms of the small parameter v = 1/T, where T 
is the relevant time scale of the problem (see Sec. HI. 
Secondly, it is an algebraic procedure that does not in- 
volve correction terms determined as solutions of differ- 
ential equations (such as Garrison's approach). Finally, 
we want a useful and practical method, one that allows 
us to do actual calculations; we want to be able to check 
the first and second order corrections formally deduced 
here against the exact solutions of many time dependent 
problems. 

To accomplish the expectations above, we need to come 
up with the right ansatz for the state |^(s)). An ideal 
ansatz should factor out the dependence of |^(s)) on all 
the terms of order 0{v^), 0(t;~^), and below. The terms 
of order 0{v~^) and below are related to e""""'-'*-' (See 
Eq. ([3])) and they are extremely oscillatory when w — > 0, 
while the zeroth order term is connected to Berry's phase 
(,i7ii(s) ^ jf w^i^ factorization could be done, we would have 
control over the divergent terms in v and immediately 
have information about the Berry phase. 

Inspired by Ponce et al. [ll| we write down the follow- 
ing ansatz for the state 1^1/ (s)) 

oo 

|*(s))=^t.P|vI/(P)(s)), (64) 

p=0 



(63) 



where 



|*(p)(s)) = ^e->"(^)e''>'"(")&(f)(s)|n(.s)) (65) 



and 



(0), 



^i^Hs) = E e^'^"'"(^)e-^"™(^)6il(s), (66) 



m=0 



with all quantities defined in Sec. [TTl We should note that 
the geometrical tcTms e'''"'''^ and e''''"'"'*' were absent in 
the original ansatz given in Ref. |1JJ. Inserting Eqs. ((65|) 
and (p^ into ([M|) we get 

oc 

l*(^))= E E^''«"""'"^'^e'^"'''^i™(«)l"(«))- (67) 

n,7n— p— 

Since the initial condition is |*(0)) = |*(°)(0)) = 



E„=o^"(0)l"(0)) it follows that 6k\0) = 5„(0) and 



|*(P)(0)) = =^ fe(f)(0) = E ^n™(0) = 0' P ^ 1- 

m=0 

(68) 
Also, imposing that the adiabatic approximation be the 
zeroth order term in the power series expansion implies 

&i°)(.) = 5(")(0) ^ b^lis) = 6^(0) = 6„(0)<5„™. (69) 



Inserting Eq. ^ into the SE, Eq. (g]), and left multi- 
plying by (fc(s)l one gets 



^ g ^Pe-^-"(^)e'^'"(^) ^A,,^is)bit{s) + b^^Js) + ijA^)b^l{s) + E A4„(.)&il(.) ) = 0. 

m=Op=0 \ 



n=0 



(70) 



Noting that 7m (s) = iAfmm(s) and 



one can rewrite Eq. (|70p in the following form 



p=0 p=0 

I 



EE^^ 



„Pe-i-^(s)e^f^(^) iA„,„(5)6(f+i)(,s) + 6(rt (,) _ M„„„(s)&il(.) + E M„,{s)b'i'i{s) 

m=Op=0 \ k=0 / 



+ ^ e->'"(^)e"^-(^)^A„™(s)6("j,(,s) = 0, (72) 



r 



where we have exchanged n <-*■ fc. The last term in 
Eq. ([7^ seems to diverge when w ^ 0. However, it does 
not because for n — m vie have A„„i(s) = while for 
n ^ m the following holds, fonm(s) = (initial conditions 
given by Eq. ((Ml) ). 

A sufficient condition to satisfy Eq. ([7^ (since its last 
term vanishes) is 

^A_(s)6(f„ti)(s) + h^^lis) + W^„™(s)6(l(s) 

+ E^^nfc(s)6i^l(s) = 0, (73) 

fc=0 

with 



W„™(s) = Mnn{s) - M„„n{s). 



(74) 



Equation ([75)1 is a main result of this paper. With the 
aid of the initial conditions given by Eqs. (p5)) and 



one can build corrections to the adiabatic approximation 
recursively. The coefficients bnm (s) are readily calcu- 
lated with the knowledge of 6„m(s), without the need 
to solve any partial differential equation as in the multi- 
variable expansion method presented in Sec. IIII A2I As 
we will show next, this fact allows us to calculate the sec- 
ond order correcting terms in a straightforward manner. 
Moreover, we have removed the highly oscillatory terms 
g--^'^r,.(s) fj-om the expression for the coefficients bnm{s)^ 
allowing a better control over any numerical algorithm 
designed to solve Eq. (j73p . i.e., it is numerically stable. 



We now proceed to calculate explicitly the first and sec- 
ond order correction terms |\l/''^-'(s)) and \^^'^\s)). The 
zeroth order term |^("'(s)) is given by Eq. (fTU)) . the adi- 
abatic approximation. 
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A. Determination of [^^^^(s)) 

When p == Eq. ^ becomes 



iA„„(s)6W (s) + 6^ (s) + M^„™(s)6W (s) 

+ J2MMs)b^Zi^) = 0. (75) 






f,(0) 



Using Eq. (jM]) we see that 6„m(s) = and that 

W„m(s)&nm(s) = W'nm (s)fo,i (0)(S„m = 0, since W„„(s) = 
0. For n ^ m the sum in Eq. ((75)l is simply Af„m(s)&m(0) 
and we get 

bl,^{s)=ih^^b„M, n^m. (76) 

When n = m Eq. (|75p is an identity and we need to 
work with the higher order expression. Setting p = 1 
and n = TO in Eq. ()73|) we have 



bl!:>{s) + J2M^,{s)b^^^{s)=0. (77) 



fe=0 



Integrating Eq. (j77p using (j76p and changing k s- m we 
obtain after using M„m(s) = — Mj^„(s), 

en^) = i^E fd^- '^^-ff 6„(0) + bW(0) 



-n-'O 



Tn=0 






(78) 



m=0 



where Eq. ([5^ was employed to arrive at the last expres- 
sion. The constant 6n„(0) is determined using Eq. 



-ir.E 



A:r,„„.(0) 

A„m(0) 



&™(o). 



(79) 



Since we now have bnm{s), for any n,™, we can insert 
Eqs. dZHl), dZni), and dBl]) into §Si) to get 

|vl/(i)(s)) = i^E e-^-"(^)e"^"(^)J™„(s)6„(0)|n(s)) 

n,7n— 
rn 7^71 

+ i^^E e-^-"(^)e'^'"(^)^=|4^'"(0)l"(^)) 

771. 7^ n 



i^E 

n,7Ti— 



e -i-.Mei7.M4M?^6„(0)|n(s)). 



Anm(O) 



(80) 



Had wc started at the GS (6„(0) = 5„o) we would get 



n=l 

+ i;,yc-^-"(^)e"'«(^)^44l"(s)) 
- i^Ee-^-^^^e-(^)^K.)), 



(81) 



where 



in the first sum. 



By looking at Eqs. ([50)1 and ([5T|l we see that they are 
identical to the ones obtained via the multi-variable ex- 
pansion method; Eqs. ([CT]) and ([55]) . respectively. Also, 
wc have a new additional term for the first order correc- 
tion, the one proportional to Jmn{s), as compared to the 
results of the standard approach, Eqs. ((22l) and (|23p . 

Using Eq. ([5T|) we can also give the conditions for the 
validity of the adiabatic approximation that comes from 
the APT by imposing that |$'^^^(s)) be negligible. 



vn'^Jno{s) 



«1, 



;SV 



'M„o(s) 



;^AA„o(s) 



, g-^i^no(s)gi7.io(s) 



MnoiOj 



A„o(0), 



<1. 



B. Determination of ^^^'(s)) 

We can proceed as before and write Eq. ((73|) for p = 1 
and n ^ 771 as 






ift 



(&i'2.(s) + M^«™(5)&i'2.(s) 



-fE^'^"fc(^)C(^))' "■^"^- (82) 



fc=0 



Using Eq. ((7^ to replace 6„m(s) above and separating 
out from the sum the term where k = m we get 



A„,„(s)V dsVA„m(s)^ 
+M„™(s)6(^)„(.^) + E M„kis)bl^l{s)).{83) 



k=0 
k^n.rn 



We can now employ Eqs. (|7^ and ([75]) to replace &nTn(s), 
n y^ m, and 6n„(s) in (|83p to finally obtain 
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6(2) (^)- (^^)^ [ d /Af„„,(s)\ ^ W,„„(s)Ai'n„(5) 
"" A„m(s) \dsVA„,„(s)/ A„,„(s) 



-Af„m(s) > — 77^Ofe(0), n 



+ Mnm(Sj 2_^ JkmXs) + 2^ ^ 7~\ f" °™(0) 



k=0 
k^m 



^mn \^ ) 



to ^"^^(0) 



^ in. 



k^n,m 



Afem(s) 



r 



(84) 



To calculate bnn{s) we set p = 2 and n = m in Eq. (|73p . 
which gives 

Cl(=^) = - E /' d.'M„„(s')6(^I(s') + 6(^2(0), (85) 

where 6Ln(s') is given by Eq. ((84|) and the constant term 
6rin(0) is determined by the initial condition in Eq. ([BS]) . 



&L'2(0) = - E ^n™(0)' 



(86) 






with &i^2i(0) obtained from Eq. ^ setting .s = 0. Fi- 
nally, the second order correction to the state |^(s)) is 

\¥^Hs))= J2 e-^-"(^)e"""(^)6(i(5)|n(5)). (87) 

n,m—0 



We should point out that, as can be seen from Eqs. 
and ([55)) . the second order correction can be calculated 
with just the knowledge of the snapshot eigenvalues 
En{s) and eigenvectors |n(s)) of the Hamiltonian H(s). 
This also holds true for the first order correction and all 
higher order terms. In other words, the APT can be seen 
as a way of converting the time dependent SE into an 
eigenvalue problem and a series expansion in the small 
parameter v. 



and 



{mr 



Mio{s)\ , T^io(s)Mio(s) 



Aio(s) ids VAio(s) 
+Mio(s)Jio(s)}. 
For a two-level system Eq. (|86)) is reduced to 



Aio(s) 



(90) 



b^ooiO) = -b^oiiO) and 6f^0) = -5^0^(0). 

Inserting the previous result into (|85p and using Eqs. 
and ([M)l we get. 



b^oJis) 



Moi(s') d /Afio(s') 



^'^^Jo ^''lAoi(s')ds' VAio(s') 



'M2 



Wwis')\A'hois ')\ 
A?o(s') 
|Mio(0)P 
A?o(0) 



'M2 



\Mio{s 

Aio(s') 



-Jio(s') 



'(i^)^ 



(91) 



and 

b?,\s) 



Mio(0)Jio(s) 



Aio(O) 
M^io(0)Afio(0) 



d /A/io(s) 



ds VAio(s) 



Aio(O) 

In the examples of Sees. I VIII and I Villi Eqs. 
be extensively used. 



s=0 

(92) 
(I92|) will 



1. Two-level system 

We now want to apply the results obtained in Eqs. 
and (|55|) to the case of a qubit. The sum in Eq. ([57)) runs 
from n,m ~ to n,m ~ 1, and the Hamiltonian H(s) is 
assumed to be non-degenerate. Thus, the second order 
correction is 

|*(2)(s)) = e-^"''(^)e'T°(")6[,o'(s)|0(s)) 

+c-^-i(^)e''^i(^'6[,i^(s)|0(s)) 

+c-i^°(«)e'''°(^)&fo^(.s)|l(.s)) 

+c-^"i(^)e'"^i(^'6fi^(s)|l(s)). (88) 

Since we assume that the qubit starts at the GS |0(0)) of 
H(0), i.e., &„(0) = 6no, Eq. dMD gives, 



^(2)r„^_ ,,^^2A^oi(s)Afio(0) 



b'i;->{s)^~iihy 



Aoi(s)Aio(0) 



(89) 



V. CORRECTIONS TO THE GEOMETRIC 
PHASE 

Let us consider a system in which its time dependent 
Hamiltonian H(s) returns to itself at the rescaled time 
Ts, i.e., H(rs) = H(0). As is well know [J], if the system 
is initially prepared in one of the eigenvectors of H(0), 
i.e, |\l/''°-'(0)) = \n{0)), and the adiabatic approximation 
is valid, then the state of the system at t^ is [^'"'(r,)) = 
ei'#'<°'(^=)|^(o)(0)). The phase 0(°)(ts) can be written as 



0(")(r3) = a(«)(r,)+/3(«)(r,), 



(93) 



where a*^'')(Ts) stands for the dynamical phase and 
P''°\ts) for the Berry phase [J: 



a(")(r,) - -c.„(s)/i;. 



(94) 
(95) 



12 



(See Eqs. (O and ([6]) for the definition of those quan- 
tities.) The Berry phase is a geometrical phase since it 
only depends on the path described by the varying pa- 
rameter in the Hamiltonian. More explicitly, if we write 
H(s) = H(r(s)), where r(s) is the parameter that is 
changed in the Hamiltonian, then [4| 7n(''"s) depends only 
on the trajectory in parameter space described by r(.s). 
For a more formal interpretation of the Berry phase in 
terms of the holonomy of a fiber bundle over the param- 
eter space see Ref. (iq . 

The concept of a geometric phase is not restricted to 
systems that start in one of the eigenvectors of H(s) or 
to adiabatic evolutions. Indeed, Aharonov and Anandan 
(AA) Q generalized the Berry phase to include those two 
possibilities. As before, we consider a non-degenerate 
Hamiltonian 13| ■ The key idea in Ref. Q was the recog- 
nition that by defining the dynamical phase as 

«(,)^_^^_ 1 rd.'(*(.')|H(.')l*(s')), (96) 
V vn In 



it is possible to show that 



/3(tc) = Htc) - aire) 



(97) 



only depends on the closed path of the curve induced by 
l^(s)) on its projective Hilbert space [13|- Here |4'(rc)) = 
ei'^(-^c)|'Jf(0)). The quantity (J){tc) is the total phase of the 
state at s = Tr and can be written as 



(rc)=Imln(*(0)|*(Tc)). 



(98) 



In the adiabatic regime, the AA phase P{tc) reduces to 
the Berry phase. Note that Tc is not necessarily the pe- 
riod of the Hamiltonian r^ . 

The AA phase is precisely the concept we need to prop- 
erly find corrections to the Berry phase in terms of the 
small parameter v defined in Sec. [l] and used in Sec. IIVI 
to build successive corrections to the adiabatic approx- 
imation. However, we need the normalized state that 
corrects the adiabatic approximation up to order p ^ j, 



with 



and 



\'${s))n,^N,\^{s))„ 



p=0 



lA^, 



(99) 



(100) 



,,, - ,(vl/(,s)|vl/(s))„ (101) 

where I^^p^s)) is defined in Eq. §5^. Following Ref. [I 
and with the aid of Eq. ([M)) we can define, up to order 
j, the following geometric phase 

/3(-'')(T,) = 0«(r,)-a(^)(T,), (102) 

where now we have 

a'-^^s) = -J^\s)/v, (103) 

0(^)(s) = lmln(*(0)|*(.s))«., (104) 



and 



fiJo 



ds\A^{s')ms')\^{s'))N,- (105) 



In our definition for /?*•■'■' we have used the period of the 
Hamiltonian t^. This is not mandatory and we could 
have chosen Tc as well. But we stick with Ts since it 
is closer to what happens in an experimental situation, 
where the Hamiltonian is slowly changed back and forth 
from its initial value. Note, however, that if r^ 7^ r^ we 
lose the meaning of /3 (Eq. (P7|) l as given by the closed 
path of l'I'(s)) on its projective Hilbert space. 



A. Zeroth order correction 

Before we show the non-trivial correction to the Berry 
phase, which is given by the first order term, it is instruc- 
tive to compute the zeroth order term. This gives us the 
flavor of what comes next without long calculations and, 
as a bonus, we are able to check that the zeroth order 
term is simply the Berry phase. We assume that the 
system starts at s = in the GS of the Hamiltonian, 

l*(0)) = |0(0)), 

although we could as well develop the same analysis for 
an arbitrary initial condition in a straightforward man- 
ner. 

The first step is the calculation of |\I'(s))7Vq, as given 
in Eq. ^ when j = 0. Since |*(o)(s)) (Eq. (HJ)) is 
normalized it is obvious that |^(s))Ar„ = |\I''-°-'(s)). Then, 
using Eqs. pUg)) and (fTU5)) we get 

a^°\^) = -^ rds'(o(s')iH(.')io(.')) 

UJo{s) 



= _J_ rds'Eois') 
On the other hand Eq. p04p gives 



(106) 



0(°) (s) = -cjo(s)/i; + 70 (s) + Im ln(0(0) |0(s)) . (107) 

Remembering that for s = r^ we have Imln(0(0)|0(T5)) = 
since |0(rs)) = |0(0)), Eq. (fTOS]) naturally leads to the 
Berry phase 



P^'\t,)= loirs). 



B. First order correction 



(108) 



We now turn our attention to the first order correction. 
As before, the first step consists in the computation of the 
explicit expression for the state |vl'(s))Arj . Using Eqs. 
and (|100[) we have 



\^is))N, = N.d^'^^s)) +v\^<''\s)) 



(109) 
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where |^'^°)(s)) and |^'^^^(s)) are given by Eqs. ((TT|) and where, without loss of gcneraUty, we have set Ni real. 
([5T|) . respectively. Had we prepared the system in an Calculating the scalar product in Eq. (|110p with the aid 
arbitrary initial state we would need Eqs. (fTO|) and ((80|) of ([81]) we get 
instead. 

To calculate the normalization constant we employ 
Eq. pm]) 



lA^iT^ ^ l + 2wRe(*(")(s)|*(i)(s))+u2(*(i)(.s)|*(i)(s)), 

where Re means the real part of a coinplex number. But 

(vl/(")(,s)|vl/(i)(,s))=in^J„o(s) 

n=l 

is purely imaginary since Jno{s) is real (Cf. Eq. (1521) '). 
Therefore, 



A^i = l/y^l + i;2(*(i)(s)|*(i)(s)). 



vi/W(5)|vi/a)(.)) = ft^(^j„o(.)) +^^E 



n=l 



?i=l 



M„o(s) 



A„o(0) 
and assuming that v is small 

7Vi = l-w2(*(i)(s)|*(i)(s))/2 + 0(w4), 



A,io(s) 
, (111) 



(110) which leads to 
I 



fll2) 



Ni = 1 



v^h^ 



T.J'^ois)] +E 



Mnois) 



A„o(s) 

r 



, g-7r'^io(s)gi7no(s) 



Af„o(0) 



A„o(0) 



+ 0(w^). 



(113) 



Notice that TVi depends on time although we have not However, the last term of Eq. (|114p can be written as 
written A'^i (s), as we have been doing with all other quan- 
tities that depend explicitly on s. Also, we have kept 
terms up to second order because they give corrections 
to first order for the dynamical phase. This can be seen 
looking at Eq. (|103p . where there exists a factor 1/v mul- 
tiplying IjJ^^'{s). 



1. Determination of the total phase 
Inserting Eq. (|109p into (|104p we get 

0(i)(s) = Imlnf(O(O)|*W(s)>+«(O(O)|^(^'(s)) 



Imhi 1 -|- iw/i ^^ J„o('''s) = arctaniift^ Jjio{ts)\ 

\ n=l / V n=l I 

n=l 

which implies that to first order 

n=\ 
= cl)(-''\Ts)+vhJ2jn0{Ts)- (115) 



where we have used Im In A^i =0 since Ni is real. When 
s = Tg we know that (0(0)|n(rs)) = Sno- Thus, 



If we use Eq. ([5^ we can rewrite the total phase as 



(0(0)|*(°)(r,,)) = e''>'«("»)e-""°("=\ 
(0(0)|*(i)(T,)) = i;ie'^°(-»)e->°(-=) Y. ^"o(^^)' 



(i)(r,) = ^(o)(r,)+«ft^/ As 



.Jo 



n=l 



A„o(s') 



116) 



and the total phase reads 



<I^^'\t,) 



Wo(t,s) 



-70 



We should note that the last term above is the first or- 
der correction to the Berry phase obtained by Garrison 
[21 and also in Ref. •3]. However, this conclusion is un- 
satisfactory for our purposes. Indeed, we are interested 
in the phase defined by Aharonov and Anandan [5|, see 
Eq. (|102p . which has a clear geometrical meaning when 
/ \ the state returns to itself (even when the adiabatic ap- 

(T^^-I-Inrln (l + ivh'S^ ^no('''.s)| • proximation fails) and is a natural generalization to the 

/ Berry phase (a]. We resolve this state of affairs in the 

(114) following. 



?i=i 



2. Determination of the geometric phase 

In order to determine the AA geometric phase we need 
to calculate, up to first order, the dynamical phase de- 
fined in Eq. (|103p . Then, subtracting it from the total 
phase computed above, we arrive at the desired AA ge- 
ometric phase. It is this first order term that we herein 
call correction to the Berry phase. 

Looking at Eq. (jlOSp we see that the first quantity we 
need to obtain is 



W (*^i)(.s)|H(s)|*(i)(s)) 



where we have used that Rc(*(°)(.s)|*(i)(s)) = 0. The 
last term of Eq. pi7p can be explicitly calculated using 
Eq. dUD, 



(vI'(i)(.)|H(5)|*(i)(.)) ^ h?E,{s){j2j^,{s) 



n=l 



h^Y.^^^') 



n=l 



Mnois) 



A„o(s) 



't'^n0(«)pi7n0(s) 



M„o(0) 



A„o(0) 
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Inserting Eq. ((TT8)) into ((TT7)) . using Eq. ((TT3)) . and keep- 
ing terms up to second order we get 



^,(*(s)|H(s)|*(s))jv, 

Mnois) 
n=l 



Eois) + 



(117) u2ft2^A„o(s 



A„o(s) 



Q-j;'^nois)^i-y„a{s) 



M„o(0) 



A„o(0) 

(119) 



(118) which, after insertion in (|103p . leads to 



a(i) (s) = a(°) (s) - «^ I] ^0 (s) - w/i- 



2 Y^ 



|Afno(0)P 



■oJno{s) + 2vh y^ Re 



A/„o(0) 
A„o(0) 



jig/g-i^„o(.')ei7.o(^') M*o(s') ) , 



r 



11201 



where we have used Eqs. ((52|) and (|106p . Notice that the 
last term has an integral of the form given by Eq. (|15p 



I 



ds'i?„o(s')e^/o''^^"^-(^"). 



with Bnoi-s) = e''^"o(*)M*o(s) and C„o(s) = -iA„o(s)/?i. 
But we have shown that this integral is at least order v 
(see Eq. (|18p ). Therefore, the overall order of this term 
is at least v^. Thus, at s = Ts, the first order correction 
to the dynamical phase is 



-^^E 



|M»o(0)P 
^1 A2o(0) 



UnoiTs). (121) 



Finally, the desired geometric phase is obtained by sub- 
tracting Eq. (I12ip from the total phase (|115p . 

/3(i)(r,) = /3(°)(T,) + 2«?i^J„o(r.) 

n=l 

+-^^E^S4lr-»o(r.). (122) 



It is worth noting that the zeroth order term above is 
the Berry phase, i.e., when w — > we have /3*-°^(ts) = 
7o(''"s) £^s our geometric phase. As mentioned before, this 
is a property any correction to the Berry phase should 
satisfy. Remembering that lj,io(ts) = 7[ /n ' A„o(s)ds and 
using the definition for J„o(s) wc can rewrite Eq. (jl22p 
as follows 

^Jo A„o(s) 



+vhY, 



\Mnom' 



f, A^o(O) Jo 



A„o(s)ds. (123) 



In Sec. IVIII we discuss how we can measure this new 
phase in general and also propose an experiment to probe 
it for the particular example of that section. 



VI. COMPARISON BETWEEN METHODS 

In previous sections we have presented four methods 
that aim to find corrections to the Berry phase and im- 
provements to the adiabatic approximation. The first 
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one, which we called standard approach, gives different 
results when compared to the multi-variable expansion 
method of Garrison [2] and the APT presented in Sec. 
IIVI However, as we have shown, to first order the last 
two methods agree. 

In the next section we show that the standard approach 
fails to properly correct the adiabatic approximation to 



first order in the small parameter v. Indeed, we show that 
the missing term in the standard approach and which 
is present in the APT is crucial if we want to have the 
right first order approximation. In other words, the APT 
developed in Sec. IIVI gives the following state for the 
time evolution of a non-degenerate time dependent sys- 
tem that starts at the GS, 



J 



l*W> 



e-^'^o(s)g 



'^°^4ll + ivh ^ J„o(.s))|0(s)) + ivhY^ 

IV n=l / n=l 



'Mnois) 



^VA„o(s) 
I 



^-i^M^)^hMs}Mjl^\n{s))} + 0{v'). (124) 



Ano(Ol 



This is the state that, to first order in v, properly cor- 
rects the adiabatic approximation. Note that it is al- 
ready normalized to first order since the normalization 
constant, Eq. (|113p . is second order in v. Furthermore, 
as we will show in the following sections, by including the 
state |^(^)(s)), as derived in Sec. IIV| we obtain the right 
second order correction. 

We have also discussed the iterative method of Berry 
[3|, who called it adiabatic renormalization [l9| because 
each iteration can be seen as a renormalization map that 
generates a new Hamiltonian from the previous one. This 
method, which is also related to other similar approaches 
[la . [l3l . Hil l , cannot be considered a perturbative correc- 
tion to the adiabatic approximation. This is because at 
each step of the iteration process v enters to all orders. 
Of course, if we stop the iteration procedure at a certain 
point we can use any method at our disposal to solve the 
transformed problem, including the APT here developed. 
In other words, we could build a hybrid approach, where 
we employ both the APT technique and the renormaliza- 
tion method of Berry. This might be an interesting topic 
to study but its full development is beyond the goal of 
this paper. 

Another method, the usual time dependent perturba- 
tion theory (TDPT), largely used to solve time depen- 
dent problems was not discussed here. The main as- 
sumption behind the TDPT is the existence of a time 
independent Hamiltonian Hq and a small time depen- 
dent part AV(i), where A <C 1. The total Hamiltonian 
is Ti{t) = Hq -f AV(i). One then builds a scries expan- 
sion in A by using the eigenvectors and eigenvalues of Hq 
(not of the snapshot 'H{t)), with the zeroth order term 
being the time independent solution to the problem. It 
is now clear what the main difference between the TDPT 
and the approaches presented in this paper is: we have 
never assumed the existence of a small time dependent 
Hamiltonian AV(t). Actually, the Hamiltonian 7i(t) can 
be seen as a particular choice of H(i), the general time 
dependent Hamiltonian used, for instance, in APT. 

We want to finish this section explaining why seem- 
ingly different approaches such as the multi-variable ex- 
pansion method and the APT of Sec. IIVI give the same 



first order correction to the adiabatic approximation. As 
we show below there is a discrete linear transformation 
that connects both approaches. This transformation can 
be written as follows 



c(f)(^,.) = J2 e-^-e'-'-^-^^bitis). (125) 

m=0 



where 7„(s) and ujjnn = ^vi ~ ^n are given by Eqs. ([5]) 
and (III), respectively Note that we will consider in the 
remaining of this section tOn as an independent variable 
{LUn{s) -^ ujn) wlicu working with expressions coming 
from Sec. IIII A2I In order to prove that Eq. (|125p con- 
nects both methods we need to show that we can go from 
Eq. dm) to dnZI) and also from Eq. (QH) to dZH) using 
Eq. (fT25l) . 

Let us start with the first part of the proof. Inserting 
Eq. p^ into ^ we get 



ri=0 p=0 



oc 



8'^'"(^)&(f^(s)|n(s)) 



Z^ z^ 

n,m—0 p—0 



„Pe->"e'^-(^)5(l(.)|n(.)>, 



which is exactly Eq. (|67|) when cum is no longer considered 
an independent variable. 

The second part requires a little more mathematical 
steps but is nevertheless as straightforward as the pre- 
vious one. Looking at Eq. (P5|) we see that it has four 
terms. We will analyze each one separately. After insert- 
ing Eq. (|125p into the first term of ([^ it results 

D^clr'-'\u;,s) ^^Y. e-^'^'""e"^"(^)A„™(s)&(f+i)(s), 

where we used £'^j(e~ """•") = — iA„„j(s)e~"'^""/?i and 
Amn(s) = — A„m(s). As is casily seen, the fourth term is 
also given by the previous expression when we set p = — 1 . 
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The second term gives the following two new terms when 
we insert Eq. (jl25p and use that ijmis) = —Mmmis), 

a,ciP)(c.,s)=^e-^-™"e'^-(^)(6(l(s)-M„„(s)6(l(s) 



Finally, after employing Eq. (|125p the third term can be 
written as 



m=0 



k,7n—0 



ijUJ^rrJlr^is) 



Mnk{.s)h^±{s). 



Putting everything back into Eq. (|29p . dividing by v, not- 
ing that iOrnn + ^n = oJm , and Considering again a;„ -^ 
iOn{s), we get exactly Eq. ^. Therefore, Eq. (fT25l) 
transforms the multi-variable expansion method into the 
APT of Sec. Ml 

Furthermore, we can also go from the APT to the 
multi-variable expansion method using the transforma- 
tion 

b^^l{s) = e-"^-(^^6„^c^P\u;,s), (126) 

where dnm ~ 1 ii n ~ m and is zero otherwise. Again the 
proof is divided into two steps. First we need to show 
that inserting Eq. (fT26l) into Eq. ^7^ we get ([27]) . 

CXD 

!*(*)) = E ^w^e-"""(")e'^'"(") 

n,m—0 p—0 

xe-'^"(^)<5„™c(f)(^,s)|n(,s)) 

OC 

n=Op=0 

which is exactly Eq. (|27p when we consider Wn(s) as an 
independent variable. To complete the proof we need to 
show that Eq. dZH), with the aid of P^ . leads to (^5)) . 
As before, we analyze separately each of the five terms 
in Eq. ((7^ . The first and the last terms are zero after 
we insert Eq. (|126p . This is the case since n = m implies 
A„„ = 0. The second term should be handled with care 
since in Eq. (jG?]) the dependence of the variables (jJ„(s) 
on s must be taken into account. This is important when 
we take the derivative with respect to s, which, according 
to the chain rule, is given by 



n=0 n=Q 



duj„ + ds 



With this in mind and remembering that — i7m(s) 
M,-am{s) we have for the second term 



E^ 

p=0 



Po-t'^"('*) 



D^c\P+^\uj,s) + dsC^^\u,s) 



where we have used the definition of Di^ given in Sec. 
IIII A2I and written out of the sum the term for p = 0. 
The third and fourth terms can easily be written as 

oo 
p=0 

and 

oo 

^„Pe-^-"(^) ^ e-^-™"(^)M„™(s)c(^)(c^,s) 

p— 771=0 

after using Eq. (|126p . In the last expression, we have also 
used that tj,„(s) = ujn{s) -\- ujrnn{s). Finally, adding all 
the terms above, multiplying the result by v, and consid- 
ering ujn{s) as an independent variable we end up with 
Eq. 



AN ANALYTICALLY SOLVABLE 
PROBLEM 



VII. 



+M™(s)c(f)(^,s))+«-V^'^"(-^)i?^cW(w,s), 



So far we have presented the general APT formalism. 
It is time to show some examples that can tell us why 
the APT of Sec. IIVI provides the right correction to the 
adiabatic approximation. For that purpose, it is desirable 
to start with a non-trivial time dependent problem that 
is exactly solved in closed form. The exact solution of 
this problem can then be expanded in terms of the small 
parameter v and compared with the results given by the 
APT. As we will see, the missing term in the standard 
approach of Sec. IIII A 11 which appears in the APT, also 
appears in the first order expansion of the exact solution. 
We also give the second order correction via the APT and 
show that it is identical to the second order expansion of 
the exact solution. We end this section comparing the 
correction to the Berry phase calculated in Sec. IVl with 
the first order expansion of the exact geometric phase 
that can be computed for this problem. As will be shown, 
both results are identical. 



A. Statement of the problem 

Let us consider a spin-1/2 (a qubit) with magnetic 
moment m subjected to a rotating classical magnetic 
field B [20]. The magnitude of the field is fixed and 
given by B = |B|. Here m = eg/{2mc)S, with e 
the electric charge of the particle, g its Landc factor, 
m its mass, c the speed of light in vacuum, and S its 
angular momentum operator. Since we have a qubit 
S = (fi,/2)cr, where cr = {iJx,(7y,(7z) are the usual Pauli 
matrices. The rotating magnetic field can be written 
as B(i) = Br{t), with unit vector written in spherical 
coordinates r(t) = {sin 9 cos (p {t), sin 9 sin (p{t), cos 9), in 
which < 9 < TT and < ip < 2tt are the polar and 
azimuthal angles, respectively. With this notation the 
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Hamiltonian describing the system is [20| 
H(i) = -m-B = 6r(t)-S, 



(127) 



where b = —Bge/{2mc) and we set e > 0. The snapshot 
eigenvectors for this problem are 

10(f)) = cos(0/2)|T)+c'^(*)sin(0/2)U), (128) 
11(f)) = sin(0/2)|T)-c''^«cos(0/2)U), (129) 

where cr^| t) = I T) and az\ I) = — | i)- The eigenvalues 
are respectively 



The important result here is that H is time independent, 
meaning that the SE for |^(f)) can be readily integrated 

|*(i)) = e""Tr^|*(0))^. Therefore, inverting Eq. (|131|) and 
remembering that |^(0)) = |^(0)), we have the solution 
to the original problem 



l*W)=c- 



■\^m- 



(135) 



Although Eq. ()135p is the general solution to the prob- 
lem, we still need to write it in a more practical way. In 
order to so, we first note that 



Eo = {h/2)b and E, = -{h/2)b. (130) fi^" = {X^ + Z^^ao and H^^+i = (X^ + Z2)"H, 



Note that the eigenvalues are time independent and we 
always have a gap of magnitude hb. 



where ctq is the identity matrix and n is a non-negative 
integer. Also, 



B. Exact solution 

If ip(t) — wt, where w > is the frequency of the rotat- 
ing ma gne tic field, the Hamiltonian (|127p can be exactly 
solved [20, l2l|. Physically, the component of the field 
projected onto the a;y-plane is rotating counter-clockwise 
around the z-axes with constant angular frequency w and 
period r = 2tt/w. This suggests that if we rotate clock- 
wise the state |^(f)), which satisfies the SE ([T]), we could 
get a new Hamiltonian H that is time independent. Let 
us define the rotated state as 



jHt 



|§(f)) = Ut(t)|*(t)), 



with 



U(f) =e~- 



lU.'t Q iwt „ 



(131) 



(132) 



where S^ = {h/2)az. Inserting Eq. (fT3T|) into the SE 
([T|) we see that |^(t)) satisfies a Schodinger-like equation 
with Hamiltonian 

H = Ut(f)H(t)U(i)-iftUt(t)^^^. (133) 

H resembles the transformed Hamiltonians of Berry's it- 
erative approach developed in Sec. IIIIBI Using Eq. (|132[) 
and the mathematical identity [20[ 

S^ cosiwt) + Sy sm{wt) = e-i^tsjhg^^in,tsjh^ 

where Sx,y = {h/2)ax,y, it is not difficult to show that 
Eq. (|133p can be written as 



H = —{bcos9 — w)az + —{bsm6)ax 

= Zgz + Xox ■ 



2 fV2 



tVH2 /iVH 



^"^ = i-UHr + U ir-Unr + --- 



4 tV4 



A^fie 



3 tl3 



5 xj5 



-l)«-(0f-(0f 



Combining both results we arrive at 

e-"^ = cos (VX^ + 2^^A ao 

, ^ sin (VX'^ + Z^l) H.(136) 

We now define three vectors that will be used later on to 
rewrite previous expressions in a more compact way. 



w = wz, 
b = fer(i), 
$7 = w b. 



(137) 
(138) 
(139) 



where z is the unity vector pointing along the z-direction. 
Since the angle between w and b is 6, the magnitude of 
Jl is simply 

^2 = w^ + |6|2 _ 2w\b\ cos e. (140) 

With this new notation Eq. (|136p can be recast as 



e R = cos I — 1 tro - ^ sm ( — ) H. (141) 



(134) 



With the aid of Eqs. (fT34)) . p35)) . plT|) . and remem- 
bering that ax\ t {[)) = I i (T))i we can calculate the 
evolution of a system that starts either at | f) or | j). 



18 






15 



— — sm f sm 



n 

nt 
T 



{b cos ( 



w) sir 




T) + 


r /m 

cos ^- 



r 



T) 
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— (6 cos — w) sin 



nt 

T 



C 2 



i), 



e^li). 



(142) 
(143) 



The most general Initial state Is written as C| | t) +cxl i), or equlvalently, 
which Implies that Its time evolution Is simply ci\4>^ {t)) + 
ci\tp^{t)). When the system starts at the GS |0(0)) of 
Eq. (|128p the time evolved state Is 



\^{t))=cos{e/2)\i'^{t)) 



-(0/2)\iA{t)), 



Mt))= 



cos 



nt 



w-h . [nt 

1 — ^ — sin — 

n \ 2 



cos(6'/2)e"'^| t) + 



e 2 



cos ( "^7 I + Q (^ cos t^ — 6) sm I — 



cos 



nt 



,w + b . /Oi 
— 1 — - — sin — 

n \ 2 



\Oit)) + '^ainesin(^]\m)}, 



sln(6'/2)e'^| i)(144) 
(145) 



r 



where, after Eqs. (|128p and (|129p . we have 

I t) = cos(0/2)|O(i)) + sln(0/2)|l(t)), (146) 

I i>=sln(0/2)e-™*|O(i)) -cos(^/2)e-™*|l(i)).(147) 

In order to avoid writing all the time |6| Instead of just &, 
we will consider & > In the rest of the paper. The final 
outcomes for all relevant quantities, nevertheless, are the 
same had we considered 6 < 0, which Is the reason why 
we win continue calling |0(s)) the GS. 



C. Expansion of the exact solution 

Since we are looking for corrections to the adlabatlc 
approximation, the frequency w = v oi the rotating mag- 
netic field should be small. An Important point Is the way 
we need to deal with terms of the form wt and w'^t. If 
we remember the definition of the rescaled time, s ~ vt, 
we see that i ex 1/u In the formalism developed for the 
APT In Sec. IIVI Therefore, the order of magnitude of, 
for example, w'^t Is the same as that of w. In general we 
have 

0(u;"+4) = 0(u;"), 

with n being an Integer. This fact should be taken Into 
account when expanding the exact solution. 
Let us write Eq. ()145|) as 

|*(t)) = no|o(i)) + ni|i(i)). (148) 

Using the definition of 51 (Eq. (|140p ) one can show that 



wcoau 



-1 + 



w sin 6 
~2^ 



0(z«3), 



which Implies that 
Ho = e"'^^ ( 1 



(^+a)t I uP' SVC? [ 



4&2 



(l_e'^«) +0{w'). 



In the previous expression, we have to expand the term 
n t. But since we now have the time t we need n up to 
third order In w 

n = b- wcose+^ sln^ 9-\ ^ cos 6* sln^ 9 + 0('w'^). 

2b 2o^ 

Using the expansion for O above and the Taylor expan- 
sion for the exponential we get 



e-^e-''"*^'"'(''/2)fl-l— sln26 

46 



,_cosesm0-^^^— 



sin'' 61 +0{W'). 



We also have the term e'^* to expand In the expression 
for Ho . But since It is multiplied by a second order term. 



„2 „:^2 



sin 0/(46^), we only need its expansion up to zeroth 



order 



gint =eibt^-iwtcos9 ^Q^^y 



Putting all the pieces together we finally obtain 



Hn 



M^_i.,.i„2(,/2) , _.U^^.^2^_^ 



46 



462 



G-{t) 



wH^ 



s\T?e + iwtcose] \ + o{w'^) 
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where 



G±(t) = l±e"'*e 



ibt^—iw t cos 6 



(149) 



Turning our attention to Hi , we see that it has an overall 
w multiplying all its other terms. Therefore, we need to 
expand 1/il up to first order 



w cost 



which results in 

Hi = e-'^-^|-^sin0(l-e'"*) 



-^sin(20)(l-c'^^*)|+O(^3). 

The second term inside the curly brackets has a nP' fac- 
tor, which means that the zeroth order expansion of e'^* 
is enough. However, the first term is multiplied by w, 
implying that we need the first order expansion of e * , 



iHi Aht—iwt cosi 



1 + 1— sm^ 



+ 0{w^). 



Using the previous expression and the expansion of 
e ' 2 up to first order we get after some algebra 

U, ^_e-i¥e--*-"'W2) /^G_(t)sin0+^sin(20) 

[_ 26 4o-^ 

X (c-it) - i^G+(t) sin 6* tan 6*) | + 0{w^). 



Finally, inserting Hq and Hi into Eq. (|148p and writing 
it as 

|^(i)) = \¥°\t))+v\¥^\t))+v^\¥^\t))+0{v^), 

we obtain 

\¥°\t)) = e-'^e-™*'^'"'(^/2)|o(i))^ (150) 

|^(i)(f)) = e-''^e-™*-"'W2)f„i!f!!!sin20|O(i)) 

\ Avb 



2vb 



G_(i)sin0|l(i)) , 



and 



(151) 



\¥^\t)) 



b t 



e 2 e 

,,2 



-iwtsin^{9/2) 



W 



4^252 



G_(i) + ^^^sin26i + iu;tcos6i) \Q{t)) 



£^sm{29)(G^it)-i^G+{t)smetane]\lit))}, 



r 



(152) 



with G±(i) given by Eq. ([MO]) . Equations (fT50ll - (fT52l) 
represent the expansions up to second order of the exact 
solution given by Eq. (|145p . 



D. First and second order corrections via the APT 

Before determining the first and second order cor- 
rections, we want to calculate explicitly the zeroth or- 
der term, namely, the adiabatic approximation given by 
Eq. pT|) . After Eq. (fTTj) one needs to evaluate two quan- 
tities: 70 (s) and uJo{s). The last one is easily obtained 
employing Eqs. (O and p30p 

ujo{s) = bs/2 = bvt/2. 

To determine 70 (s) we need Mqq{s) as given by Eq. ([5]). 
Using Eq. (|128p for the snapshot eigenvector |0(s)) we 
get 



|6(s))=i^sin(0/2)e'^U), 



which implies 



Moo{s)^i-sm\e/2). 

V 



(153) 
(154) 



Thus, inserting Eq. ()154p into ([5]) we get 



7o(s) = ssm\e/2) 

V 

and Eq. pT|) reads 

|vl/(«)(i)) 



-wt sm\0/2), 



e 2 e 



■iiutsin^(6l/2) 



(155) 



\0{t)). (156) 



The first order correction obtained via the APT in Sec. 
HVl is given by Eq. ([8T|) . Since we deal with a two-level 
system, there is no sum and we can set n = 1 in all terms 
of Eq. ((8T|) . In addition to 70 (s) and ujo{s), we need 
to compute 71 fs), wi(s), A/io(s), Aio(s), and Jio(s) to 
determine |^'^'(.s)). Wc start with the gap, which is 
easily computed using Eq. p30p 

Aio(s) = -hb. (157) 

Using Eqs. (O and (|130p we immediately get 

LUi{s) ~ ~ujq(s) ~ ~bvt/2. 

The term 71 (s) is obtained after calculating Mii(s). Us- 
ing Eq. (fTJgi) we get 



|i(^)) 



-i-cos(0/2)c'-U). 

V 



(158) 
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Inserting Eq. p58)) into ^ wc find that 
Mii(s) =i-cos2(6'/2), 

V 



(159) 



which leads to 

71 (s) = -wtcoa'^{e/2). (160) 

Using Eqs. ©, P^ . and p55)) wc arrive at 

Mio(s)--i— sin0. (161) 

2v 



FinaUy, with the aid of Eqs. ^, (fTST]) . and plT|) we 
obtain 



Jiois) 



Avhh 



srn^e. 



Therefore, returning to Eq. ((8T|) using that 



Afio(s) . w . , 
~ i „ .^ smt 



Aio(s) 2i;6n 



(162) 



(163) 



we get 



|vl/(i)(t)) = e-''^e-™*^'"'(''/2)('_i^sin20|O(i)) 



2w6 



G_(f)sin0|l(t)) , 



(164) 



where G-{t) is given by Eq. (|149p . 

Moving on to the second order term, Eq. ((88|) . whose 
coefficients are obtained from Eqs. ((89|) to ([92|) . we 
sec that almost everything we need to explicitly write 
\^^'^\t)) is already calculated. We are left with only two 
quantities to compute, which are 



d / Mio(5) 
ds V^io(s) 



0, 



as can be seen from Eq. (|163p . and 



w 



Ww{s) = i— cos6' 



(165) 



where we have employed Eqs. ([741), (|154[) . and (|159|) . We 
are now able to write down explicitly the values of the 
four coefficients. The first one, Eq. ([89|) . is easily calcu- 
lated by noting that Afoi(s) = Afoi(O) == -Mj*'o(0). Thus, 






(166) 



The second one, Eq. ([90]) , is obtained inserting the values 
of W^io(s), Mio(s), Aio(s), and Jio(s), 



w^sin(26') / wt 



&ro^(^) 



1 — i^— sin Stan! 



(167) 



The evaluation of the third coefficient, Eq. (PT|) . is just 
a little more involved. The integrations are easily done 
since the first integrand is time independent and the sec- 
ond one is a linear polynomial of the rescaled time s. 
Putting the results of the integration back into Eq. (pij) 
we can rearrange it as follows. 



.f2i/ N w sm\ , 



w'^t^ SlT?t 



-\-iwt cos ( 



(168) 
The fourth and last coefficient, Eq. ((5^ . is calculated 
in the same manner as we did for fejQ (s). After some 
algebra we get 



,f2)/ X w^sin(26') /, .wt . „ 
Ai s = , ,,, ' 1 + i-- sin 61 tan ( 



(169) 



Inserting all the coefficients above into Eq. ([881) we get 
after some algebraic manipulations 



\¥^){t)) 



b t 



e 2 e 

,,2 



-\wts\n^{e/2) 



W . 2/ 

■ sm 6 
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G_(i) + ^^^sin26' + iu;tcos6') |0(f)) 



^ sin(20) f G_(t)-i^G+(t) sine tan 0)|l(t))^, 



r 



(170) 



with G±(t) given by Eq. ^^. 

Wc are now in position to reach interesting and im- 
portant conclusions. First of all. comparing Eqs. (|150p - 
(fT52l) with Eqs. (fTSe]) . (fTei]) . and (fTTO]) wc easily reahze 
that they are the same. In other words, the expansion of 
the exact solution up to second order is identical to the 
correction to the adiabatic approximation up to second 
order obtained from the APT of Sec. IIVI Second, since 
Eq. psip and (|164p agree, we can rule out the standard 



approach of Sec. (jIII A ip as the right way of correcting 
the adiabatic approximation. Indeed, the term propor- 
tional to 

-i^sin20|O(i)) 

is absent in the standard approach first order correction. 
Although not shown here, we also obtain different second 
order terms whether we use the standard approach or the 
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APT. And evidently, the correct term comes from the 
APT, as Eqs. p^ and pTD)l demonstrate. 



E. The geometric phase 

We have demonstrated in the previous paragraphs that 
the APT gives the right first and second order correction 
terms to the adiabatic approximation. In this section our 
goal is to prove that the formalism developed in Sec. |Vl 
and which rests on the APT, is also the appropriate one 
when one is interested in corrections to the Berry phase. 
We first need to calculate the exact geometric phase for 
the state given by Eq. ()145p . We then expand this phase 
in terms of the small parameter v = w, allowing us to 
compare it with the first order correction obtained via 
the formalism of Sec. IVl 



1. The exact geometric phase 

We are interested in the geometric phase that the state 
in Eq. (|145|) acquires after the Hamiltonian H(i) returns 
to itself. Looking at Eq. (|127p we see that the period of 
the Hamiltonian is r = 2tt/w, or r, = 2'kv/'w if we work 
with the rescaled time. The geometric phase we want to 
calculate is given by Eq. (P7|) . Therefore, we need first 
the total phase (J){t) and the dynamical phase a{T). 

The total phase, Eq. (|98)) . is obtained using Eq. (|145p . 
which gives the state of the system at t = t. At t ~ 
0, on the other hand, we have |^(0)) = |0(0)). Hence, 
remembering that {n{0)\m(T)) ~ Snm we get 



(*(0)|*(r)) = 


^ p-™r/2 


cos 

V 2 , 






. w cos 9 — b . 


v2 j. 


; 


= g-i»r/2^ 


? e<, 





with R = |(*(0)|*(r))| and C = arctan (lm(*(0)|*(r))/ 
Re(^(0)|*(T))). Therefore, using Eq. dMl) we get for the 
total phase 0(t) = —wt/2 + (, or more explicitly 



WT 



arctan 



iO-b 



n 



tan 



~2 



(171) 



The dynamical phase is given by Eq. (j96p . which in 
terms of t is 

a{T) = -\j\\t{-^{t)mt)\^{t)). 
Using the definition of il, Eq. ()140p . we get 



hb 



2w' 



^{t)\n{i)\^{t)) = - 1 - --sin^0sin^(m/2) 



^2 



which results in 

br w^bTsin^O 



a{T) 



u)26sin(J7r)sin^6' 
2?F ■ 



(172) 



The exact geometric phase, Eq. ((97)) . is calculated sub- 
tracting from the total phase the dynamical phase. Thus, 
using Eqs. pTTj) and (fTTS]) wc get 



/3(r) 



WT 



arctan 



br 

y 



w'^brsin^O 
2f^2 



w cos 9 — b ( fir 

— n — *"H~ 

w'^bsin{nT)sm'^6 



2f}3 



• (173) 



2. Expansion of the exact geometric phase 

We now proceed with the expansion of the exact results 
obtained above up to first order in the small parameter 
V = w. Again, we should be careful when doing such an 
expansion since we are always assuming to be near the 
adiabatic regime. This implies that the period r of the 
Hamiltonian is a large number of order l/w. Therefore, 
terms like w'^t are actually 0{w), which means that we 
need to expand all expressions up to second order in w 
and then look after terms of this type. 

Let us begin with the total phase. Using the definition 
of 51 and expanding the inverse of the tangent given in 
Eq. (J17ip wc obtain up to second order in w, 



c 



br wTCOsB w^T sin 6 ly^ sin 9sm{bT) 
Y ^ 2 46 ^ 462 



The last term is second order in v since | sin(6r)| < 1, even 
for large r. The other term containing w'^ is, nevertheless, 
0{w) because it is multiplied by r. Hence, the total phase 
expanded up to first order is 



, , br . ,,^ ,„, ii;2Tsin i 

(t) = WT sin^(6'/2) - 



46 



Oiw^). (174) 



The dynamical phase up to first order is obtained not- 
ing that the last term of Eq. (|172p is 0{w'^) since 

9 9 

^ sin(17r) = ^ sin(6r) -f ©(u.^). 
il"^ b'^ 

Then, using that w^/fl'^ = w'^/b'^ + 0{w^) wc get 
6t w'^t sin '' 



a(r) 



26 



+ 0(u;2), 



(175) 



which leads to the first order expansion of the geometric 
phase below, 

(3{t) = -WT sm\e/2) - ^'"'^f'''^ + 0{w^). (176) 

46 



3. Perturbative correction to the geometric phase 

As shown in Sec. |Vl the zeroth order term of the ge- 
ometric phase defined in Eq. p02p is simply the Berry 
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phase. For the particular problem of this section it can 
be easily calculated using Eqs. ([5]) and (|154p . 

/3(o)(t,) = -WT sin\e/2) = -wt{1 - cos6')/2, (177) 



where we have used that r, = vt 
T we get /3*^°^(ri5) = — 7r(l — cos 6*^ 



Using the value for 
This phase can be 
interpreted as half of the solid angle subtended by a curve 
traced on a sphere by the direction of the magnetic field 
while it goes back and forth to its initial value Q . 

The first order correction to the Berry phase is calcu- 
lated by using directly Eq. p22p , 



Inserting Eqs. (|162p . (|163p . and noting that loiq{ts) = 
— 6t, = —bvT we get 



l3'^^\Ts) = -WTsm\e/2) 



3wVsin^6l 
4& ■ 



(178) 



Comparing Eq. (|178p with the expansion of the exact ge- 
ometric phase given in Eq. (|176p we see that they are 
identical. In other words, the previous result shows that 
we get the same answer for the correction to the Berry 
phase either if we expand the exact AA geometric phase 
or if we calculate the AA geometric phase for the correc- 
tion to the adiabatic approximation given by the APT. 
However, and it is here that the usefulness of a pertur- 
bative method becomes evident, for the vast majority of 
problems we do not know their exact geometric phases 
and we must rely, therefore, on the APT and the methods 
of Sec. |V]to go beyond the Berry phase. 



^. Measuring P^^' {ts) 

The correction to the Berry phase j3^^\ts) can be mea- 
sured as follows. We prepare a beam of particles in the 
GS |0(0)) of the Hamiltonian H(0) and split it into two 
equal parts. Half of it is subjected to the time dependent 
Hamiltonian H(s) and the other half to a time indepen- 
dent one, H(0). In the first beam H(s) is changed with 
time in a manner that makes the first order correction 
to the adiabatic approximation relevant. This is done by 
adjusting the frequency w of the rotating field. For the 
other beam, H(0) is such that it gives the state |$(ts))jvj 
= e'"('^=)|\I'(rs))Ari at s = r^, i.e., the state |$(rs))Ari as 
given by Eq. (|109p with an additional phase equals to 
the dynamical phase of |\I'(rs))jVi. This is achieved by 
|^('''s))Afi being an eigenvector of H(0) with an eigen- 
value set in a manner that provides the phase a(rs) at 
Tg. Then, recombining the two beams we measure its in- 
tensity for several orientations of the magnetic field (the 
angle 9). An interference pattern emerges whose inten- 
sity contrast is proportional to cos^(/3'"'^-'(6')), which can 
be compared with the contrast predicted by Eq. (|178p . 



It is worth noticing that it may not be easy to build ex- 
perimentally the Hamiltonian H(0). 

We want to end this section analyzing the case where 
Tg = Tc, i.e., where the periodicity of the Hamiltonian 
Ts = 2ttvIw equals the time that it takes for the initial 
state to return to itself up to an overall phase [a]. In Sec. 
IVl we emphasized that those two periods are in general 
different. If one looks at Eq. (|145p it is straightforward 
to see that the exact solution returns to itself (up to an 
overall phase) after a time Tc = vt = 2TTv/il,. However, 
in general we do not know the exact solution and we must 
rely on the period for the corrected state to return to itself. 
To first order the system comes back to the initial state 
when the term multiplying the state |l(s)) is zero at s = 
Tg — VT. From Eq. (1164^ this is the case when G_ (r) = 0, 



-Bge 



l + cos6' 2mc(l + cos( 



(179) 



after using the values for r and b. Since w <^ 1 this condi- 
tion can be achieved by choosing a small field. If possible, 
we can also choose a particle with either a small charge 
or a big mass, or change the orientation of the field. But 
assuming this condition is fulfilled the geometric phases 
defined in Sec. |V] acquire the geometrical meaning that 
is inherent to the AA geometric phase Q. 

Indeed, using Eq. (|179p and r = 2Tr/w, the first order 
correction to the Berry phase given by Eq. (|178p becomes. 



/3(i)(r,) = -27rsin^(e'/2)- 



2(1 -h cos 6*)' 



(180) 



which only depends on the angle 9, i.e., the angle of 
the magnetic field with the z-axis (there is no other dy- 
namical component here such as the small parameter 
V ~ w). Employing Berry's phase definition we can write 
Eq. HIDD as 

/3(i)(r,) = 7o(t.) + ^7o(r.) = ^7o(r.). (181) 

This is the geometric phase when the first order correc- 
tion to the adiabatic approximation is relevant, and it 
can be probed by using an experimental setup similar 
to the one developed to test Berry's phase Q with the 
following slight modification. 

First, a polarized beam of spin-1/2 particles prepared 
in the GS |0(0)) is split into two beams that are sent to 
regions with magnetic fields pointing initially in the same 
direction (see Fig. ^. In one path the direction of the 
magnetic field is kept constant and its magnitude {B2) 
is tuned such that at s = r^ the phase of the particles is 
given by the dynamical phase a*-^-* (ts) (Eq. (|175p ) with w 
satisfying Eq. (|179p . Note that the GS is independent of 
the field strength (Eq. (|128p ). Along the other beam the 
field (i?i) is slowly rotated with frequency uj back and 
forth around the z-axes. The frequency should be con- 
sistent with p79p and chosen in a way that makes the 
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FIG. 2: (Color online) A beam of particles prepared in the GS 
is split into two equal parts. One (upper beam) goes through 
a region of constant magnetic field whose strength B2 is such 
that at the end it acquires the dynamical phase a^^' (tb) of the 
lower beam. The latter beam goes through a region where the 
magnetic field Bi rotates around the z-axes until it returns to 
itself. (In the original proposal [J| , the field strengths are the 
same, Bi — B2-) Finally, the beams are recombined and the 
intensity measured, allowing us to determine the geometric 
phase /3'^'(rs). See text for more details. 



first order correction to tlie adiabatic approximation rel- 
evant. Then the beams are recombined and the intensity 
measured. Repeating this experiment for several values 
of 9 we should sec the intensity changing as cos^(/3(^^ (6^)), 
where (3^^\6) should agree with Eq. (J18ip . 



VIII. NUMERICAL EXAMPLES 

In this section we want to consider three more exam- 
ples and compare their exact time evolution with the first 
and second order corrections to the adiabatic approxima- 
tion given by the APT. One of the examples can be seen 
as a particular case of the analytic problem in Sec. IVIII 
and another one can also be solved analytically in terms 
of a special function (See Appendix E]) . However, here 
we solve them all numerically. 

We again restrict ourselves to a two-level system de- 
scribed by the following Hamiltonian, 



H,(s) = 



£:c'''^("^) 



(182) 



where 2E is the time independent gap of the system and 
Oj{s) is the time dependent part of the Hamiltonian. We 
choose three polynomials for Oj{s), j = 1, 2, 3, which de- 
fine our examples: 



9,i.s) = 9^ + w,s- 



(183) 



The parameter 9'^ represents the initial condition for 
9j(s) and Wj > 0. For j = 1 we recover the example 
of Sec. IVIII when the angle of the magnetic field with the 
2;-axes is 7r/2. Note that we are already working with the 
rescaled time [22[ . 

The snapshot eigenvectors and eigenvalues of Hj (s) are 

|0(s)) = 4= fe'^^(^)| t) + I i)) with Eo = E, (184) 



V2^' 



|1(5)) 



1 

71 



U9,{s)\ 



|t)-U) with £:i = -^(185) 



An arbitrary state at s can be represented as 

l*(s)>=CT(.s)|t)+q(s)U). (186) 

with coefficients satisfying (e = E/{vh)) 

ct(s) = -i6e"'^-(^)cx(s), (187) 

cx(s) = -iee-"'^(^)cT(s). (188) 

The comparison between the exact time evolution of 
l'I'(s)) and the approximate results of the APT simplifies 
if we rewrite Eq. (|186p in terms of the snapshot eigenvec- 
tors of H(s). Using Eqs. PM)) and (fT^ onc gets 

It) = e-«^(^)(|0(.)) + |l(.)))/V2, (189) 

li) = (|0(.))-|l(.)))/^/2, (190) 

so that (|186p becomes 

\^{s))=co{sms))+c,{s)\l{s)), (191) 

where 

coi-s) = (e-'«^■(^'cT(s)+q(s))/^/2, (192) 

ci(s) = (c-'«^(^)cT(s)-q(s))/V2. (193) 

If the system starts at the eigenvector |0(0)), i.e., co(0) = 
1 and ci(0) = 0, then 



c? 



(0) = e"'?/V2 and C|(0) = l/%/2 



(194) 



To have a quantitative measure of the closeness of the 
corrections to the adiabatic approximation to the exact 
state (|19ip we compute a quantity called fidelity. 



F,{s)^\{^{s)\^{s))nJ', 



(195) 



where |\I'(s))Arj, is the normalized state containing cor- 
rections up to order k (Eq. ([99]) ). When the states are 
the same Fk ~ I and Fk = when they are orthogonal. 
Using the snapshot eigenvectors given by Eqs. p84p 
and (|185p . Eqs. 199]) and (|10ip . and repeating the same 
steps of Sec. IVIII we get 

|*(s))jv„ = |*^°Hs)> = e-'"e-''^«^(^)/2|0(s)), (196) 
with A9j{s) ^9j{s)-9'^j, 

\^{s))m = Ni (l^^'Hs)) + z;|vI/W(5))) , (197) 

where 



e-ieSg-iA0,(s)/2 J _ 



ih 



8E 



9j{s')ds'\0{s)) 



+ — (^,(s)-e'2-0,(O))|l(.))[, (198) 



and 



(199) 
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in which |$(^)(s)) is given by Eq. (|88|) . The coefficients 
of [^^(^'(s)), where wo(s) = — Wi(s) = E s/h and 7o(s) = 
71 (s) = -A6'j(s)/2, are 



b'^K^) = "'' ^'^^ 



^^^ iQ) + e]{s) + -^[l e]{s'W^ 



C^(5) 



&fo^(.) = -'^ 



0,(O)0,(s), 



^^ ^ (.) ^ ^^(^) 




4 / "?^"'^^"' 



]{s')ds' 



^^hs)-^A-om-'-^ 



e|(s')ds' 



By inspection of Eqs. dUS]), (fTOTl) . P^ . and their coef- 
ficients, and using the definition for Oj{s), we reahze that 
from one order to the next we have a smaller contribu- 
tion to the overall state if e""'^ = vh/E < 1. The previous 
condition is related to the existence of a gap {E > 0) 
and the near adiabaticity approximation {v = Wj <^ 1). 
When those conditions are satisfied, we should expect the 
APT to work. 

There is one more interesting fact. If we factor out the 
highly oscillatory dynamical term e^"^**, the other oscil- 
latory terms are always multiplied by the first or second 
order derivatives of ^j at s = 0. This can be seen by 
looking at Eq. (|198p . where we have the term &^'^''9j{0). 
A similar exponential appears in |^'^)(s)), multiplying 

(2) (2) 

either 6*^(0) or 6*^(0) (see coefficients 6oi (•^) ^^'^ ^ii (•*))• 
Therefore, by properly choosing the functional form of 9j 
we can eliminate those oscillatory terms. It remains only 
a global oscillatory phase e""^** that has no influence on 
the fidelity or on the probability to find the system out 
of the GS. 

Let us start presenting the results of the numerical 
calculations. In Fig. [3] we show the value of the infidelity, 
|1 — Fk{s)\, when e^^ < 1. For the three cases, as we 
increase the order of the APT we get closer and closer to 
the exact solution (small infidelity). In Fig.[3]we show the 
behavior of the APT as we increase e^^. We computed 
how much the second order correction differs from the 
exact solution for all Oj{s). It is clear that for e^^ < 1 
we almost sec no difference from the exact solution. For 
e~^ > 1, however, the perturbation theory fails as can be 
seen from the last panel of Fig. 21 

In all previous calculations it was implicit that 0j{s) 
was a smooth function. It may happen that its first or 
second (or n-th) order derivative with respect to time 
becomes discontinuous. This is related to the way we 
can experimentally control the Hamiltonian [llj. Under 
those circumstances we can continue using APT to pre- 
dict the behavior of the exact solution to the SE. The 
way to circumvent this problem is relatively simple. Let 
us assume we have the following functional form for 0j{s) 
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FIG. 3: (Color online) Here 9° = 1, £■ = 2, and u = Wj = 0.5, 
which gives e~^ = 0.25 (ft = 1). At the top we have 9j{s), 
j — 1,2, and at the bottom j = 3. The black/dotted curves 
represent the infidelity between the zeroth order correction, 
Eq. (|196|l . and the exact solution, Eq. (|19ip . as a function of 
the rescaled time s. Both quantities are adimensional. The 
blue/dashed curves are the infidelity when we go up to first 
order (Eq. (|197|l ') and the red/solid ones when we include the 
second order term (Eq. (|199p ). For j — 1, the first and second 
order curves are indistinguishable and the solid/dotted curves 
go as high as 0.004. 
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FIG. 4: (Color online) The same parameters of Fig. |3] but 
with difi'erent gaps. Top: e"^ = 0.125 and 0.25. Bottom: 0.5 
and 5. All curves represent the infidelity between the exact 
solution and the adiabatic approximation corrected up to sec- 
ond order (Eq. (|199|l V The solid curve represents 9-i{s), the 
dashed ^2(5), and the dotted one 9i{s). In the first panel all 
curves coincide while at the next two the dashed and dotted 
curves are indistinguishable. Note the difference of scale at 
the bottom panels. For the first three, the APT works beau- 
tifully and the results are better the lower e^^. At the last 
panel we see the three curves and the break down of the APT 
since e~^ > 1. 



When s < 0, and starting, let us say, at s = —0.2, 
and using the initial condition at that time, we com- 
pute the perturbative terms given by the APT using 
9j{s) = 6*°. All terms but the zeroth order vanish since 
the Hamiltonian is time independent for s < 0. Then, at 
s =: we start computing the perturbative terms using 



,(^) 



90 



WjS^ and as initial state we use the final 



25 



state from the previous computation, i.e., wc impose the 
continuity of the wave function at s = 0: Um^_,o- l^ls)) 
= lim5^o+ l*(-5))- This procedure allows us to obtain in 
a perturbative way the right time evolution for the whole 
range of rescaled time s. We exemplify this approach in 
Fig. [5l It is clear that this approach (third panel) is the 
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FIG. 5: (Color online) The same parameters and notation 
of Fig. |3]for the case 02 (s). At the top panels we used, for 
expressions coming from the APT and throughout the whole 
range of s, 0(0) = for the first panel and 0(0) = 2w2 for 
the second one. At the bottom panel we used 0(0) = for 
s < and 0(0) = Iw^ for s > plus the continuity of the 
wave function at s = 0. 

best option. To second order, we see no appreciable dif- 
ference between the exact solution and the perturbative 
solution. In Fig. [5] we worked with O^i^s) but the same fea- 
ture shows up with 0i(s), where in this case it is the first 
order correction that is problematic. The same feature 
is true if wc work with another time dependent Hamil- 
tonian. In general, a discontinuous derivative of order 
k + 1 in the quantity Mnm{s) I ^mn{s) affects the k-th or- 
der in perturbation theory. The remedy, nevertheless, is 
the same as before. 



IX. DISCUSSIONS AND CONCLUSIONS 

In this paper wc presented a useful and practical way to 
find corrections to the adiabatic approximation named, 
after Garrison [31, adiabatic perturbation theory (APT). 
Considering the adiabatic approximation as the zeroth 
order term, wc have developed a power series expansion 
that gives the time evolution of the system. The only 
assumption made was the existence of a non-degenerate 
Hamiltonian throughout the time evolution. We have ex- 
plicitly calculated corrections up to second order in the 
small parameter v, that is related to the inverse of the 
relevant time scale of the problem, namely the time re- 
quired to change the system's Hamiltonian from its initial 
value to the desired final one. 

We have checked the validity of this approach compar- 
ing the exact solution of several time dependent prob- 
lems with the approximate results given by the APT. One 



of the problems had an exact analytical solution which 
allowed detailed comparison with the approximate one 
given by the APT. We got a perfect agreement between 
both ways of solving the problem. The other time depen- 
dent problems were solved numerically. The APT passed 
all tests for those numerical cases too: the more terms 
one adds to the approximate solution the closer one gets 
to the exact solution. We should note, however, that a 
rigorous general proof of convergence of the APT series 
expansion was not given, although we believe that it will 
work in general at least in an asymptotic sense. 

In addition, we have compared the APT to other meth- 
ods that also try to go beyond the adiabatic approxima- 
tion. The first method we dealt with was what we called 
the standard approach, since it is based on the straight- 
forward manipulations of the integral equations that one 
gets when writing formally the exact solution to the time 
dependent Schrodinger equation (Sec. IIII A ip . Wc have 
shown that the naive expansion of the integral equations 
in terms of the small parameter v fails to give an accurate 
correction to the adiabatic approximation. We then stud- 
ied the iterative rotating-basis method developed in Ref. 
3 and which is related to the ones in Refs. |12. [l3l. Il4| . 
As can be seen in the analysis of Sec. IIIIBl this approach 
is not a perturbative method in the small parameter v. 
Rather, it is built on another premise that, loosely speak- 
ing, has the goal of finding by an iterative process a new 
frame of reference where the modified Hamiltonian be- 
comes time independent. We have emphasized that at 
each iteration step one can in principle use our APT as 
a way of approximating the solution within that frame. 

Most importantly, we have proven that the APT here 
introduced, and which was inspired by the work of Ponce 
et al. jll| . is connected to the multi- variable expan- 
sion method developed by Garrison 01 . Indeed, we 
have shown the formal mathematical equivalence be- 
tween both methods. Starting with the APT wc can ob- 
tain the multi-variable expansion method and vice versa. 
However, the equations obtained from the APT to or- 
der p are simple algebraic recursive relations involving 
the terms of order p — 1. On the other hand, the multi- 
variable expansion method requires not only manipulat- 
ing recursive relations but also solving partial differential 
equations. 

We have also shown how to calculate corrections to the 
Berry phase [^ to an arbitrary order in the small param- 
eter V. The strategy we adopted had two basic ingredi- 
ents, one of which was the normalized p-th order correc- 
tion to the adiabatic approximation. The other one was 
the Aharonov-Anandan phase, a natural generalization of 
the Berry phase |5j , suited to the calculation of geomet- 
ric phases away from the adiabatic regime. Moreover, we 
have explicitly computed the first order correction in a 
spin-1/2 (qubit) problem, and proposed a specific quan- 
tum interference experiment to measure it. We showed 
that when the first order correction to the adiabatic ap- 
proximation is relevant, the geometric phase should be 
two and a half times the Berry phase. 
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Finally, our results lead naturally to new questions. 
First, can we build an APT similar in spirit to the one 
presented here but for open quantum systems where we 
have non- unitary dynamics [23|? Second, can we employ 
this open dynamics APT to calculate corrections to all 
sorts of geometric phases 12J1? And third, can we extend 
our ideas to the case where the Hamiltonian spectrum is 
degenerate? 
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APPENDIX A: SOLUTION TO THE 62(3) CASE 



For 6*2(5) = 0°2+W2 s^ Eqs. p87)) and p88)) are a par- 
ticular case of the following ones, 

c^{s) = -F*e-™-'cT(s). 



Decoupling we get, 

ct(s) - i2w2sc^{s) ~ \V\'^c^{s) = 0. 

Making the change of variable c^{s) ~ f{s)z{s) 
and imposing that in the new equation the coef- 
ficient multiplying z(s) be zero we obtain z(s) -t- 
{iw2+wls^ + \V\^) z{s) = 0, with f{s) = /(0)e™^'*'/2^ 
Making another change of variable, x — ■\/2|w^s, we get 
(Pz/dx"^ + (a;V4- a) z{x) = 0, where a = -\V\^ - i/2. 
The solution to the previous equation are the Weber func- 
tions (ii, 



zi(.t) 

Z2ix) 



a2„:r2"/(2n)!, 



n=0 

00 

n=0 



7(2n+l)!, 



in which oq = ai = 1, 02 = 03 = a, and an+2 ~ cicin — 
n(n— l)a„_2/4. Finally, returning to the original variable 
we get the solution to the original problem. 

c^{s) = e™^^'/2 (cOzi(2/M s) + c°Z2(2/M s) 

with c" and C2 being fixed by the initial conditions. 
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